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ABSTRACT 

The  present  study  was  an  experimental  investigation  of  the 
surface  heat  transfer  distribution  about  cones.  The  study  was  conducted 
in  the  Princeton  University  hypersonic  nitrogen  facility  covering  a 
wide  range  of  geometries  and  test  conditions  in  an  attempt  to  evaluate 
the  numerical  accuracy,  theoretical  modeling,  and  range  of  applicability 
of  the  several  theories  which  are  available.   Cones,  varying  from  3° 
to  30° )   with  two  nose  bluntnesses,  and  lengths  from  one  to  eight  inches 
were  examined  over  a  range  of  test  conditions,  Mach  numbers  from  16  to 
25  and  unit  Reynolds  numbers  from  7,000  to  40,000  per  inch.   The 
detailed  heat  transfer  distributions,  covering  the  full  range  of  "nose 
dominated",  "merged",  and  "weak  interaction"  regimes,  were  compared 
with  several  theoretical  prediction  techniques.  A  critical  review  of 
the  modeling  involved  in  the  theories,  as  compared  to  the  experiments, 
is  presented. 

The  investigation  indicates  that  the  effects  of  transverse 
curvature  and  bluntness  have  been  overestimated  in  the  past.  The 
classical  boundary  layer  theory  is  shown  to  predict  surface  heat 
transfer  quite  satisfactorily  over  most  of  the  range  of  the  study. 
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v 

A  5   cos  9 

5  Boundary  layer  thickness 
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Blunt  cone  coordinate  defined  by 
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Flow  deflection  angle 


vi 


9  Cone   half  angle 

c 

9  Shock  angle 

s 


[i  Viscosity 

p  Density- 


CD 


Shearing  stress     =     \i  (— ) 

dy  y  =  o 

Hypersonic  boundary  layer  interaction  parameter 

-IET 


=  M3    W, 


Re. 


\jr  Stream  function  such  that  L      =  2jtrpu 


^r  Non-dimensional  form  of  stream  function 


* 


2 

p  uitr 

00     00 


a)   ,ai.  Constants   in  tangent  cone  formula 


=  7(7+1X7+7) 

(7+3)' 


,2 


2 
3  7     -  77+3 

^1  "  2(7+3)(7+l) 


CD 
1  1 


1  U)     ,,  2,      2 

o  M     tan  0 

oo  c 

Subscripts 

w,oo,o         pertaining  to  wall,  freestream,  and  stagnation  conditions 
respectively 

b  refers  to  spherically  blunted  cone 

c  pertaining  to  inviscid  cone  flow 

e  pertaining  to  conditions  at  the  edge  of  the  boundary  layer 

M  denotes  classical  Blausius -Mangier  solution 
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s  refers  to  pointed  cone 

r  Recovery  condition 

T  Pertaining  to  nozzle  throat 

*  Based  on  reference  temperature  T 
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INTRODUCTION 

The  present  study   is  an  experimental   investigation  of  the   heat 
transfer  distribution  about  cones.     With  the  advent  of  ultra-high 
speed  flight,   an  ability  to  predict  local  heat  transfer  distributions 
in  rarefied  flows  has  become   increasingly  more   important  to  the  design 
engineer.      In  a  broader  sense,   the  present  work  is  part  of  the   con- 
tinuing program  of  theoretical  and  experimental  investigation  at  the 
Gas  Dynamics  Laboratory  which  has  been  focused  on  the  problem  of  the 
interaction  between  the  boundary  layer  on  a  slender  body  and  the 
surrounding  hypersonic   stream.     Although  the  measurement  of  surface 
parameters   such  as  heat  transfer,   pressure  and   skin  friction  are 
insufficient  in  themselves  to  provide  a  thorough  understanding  of  the 
problem,   it  is  intended  that  the  present  heat  transfer   investigation 
will  reinforce  other  studies  which  have  concentrated  on  the  displacement 
and  structure  of  the   shock,   the  viscous  layer  profiles  and  other  details 
of  the  flow-field  development. 

The  experimental  heat  transfer  results  were  obtained  for  slender 
cones  whose  half -angles  were  3°  ,   5°  ,   10°  ,  15°  ,   20°    and  30°  .      In 
addition  to  the  extensive  sharp  cone  studies,  blunted  cones  of  the 
same  half -angles  (excepting  the  3°    cone)  were  investigated  to  determine 
the  effect  of  the  nose  region  on  the  conical  afterbody.     Tests  were 
conducted  over  a  Mach  number  range  from  16-25  and  unit  Reynolds  numbers 
from  7,000  to  40,000  per  inch. 

The  flow  field  over  a  cone  encompasses   several  fluid  mechanical 
flow  regimes  ranging  from  the  classical  thin  boundary  layer  and 
inviscid  outer  flow  at  large  distances  downstream  to  the  free  molecular 
flow  at  the  vertex.     To  evaluate  the  numerical  accuracy  and  determine 
the  validity  of  the  various  theoretical  flow  models  proposed,   the 
results  are  compared  to  selected  theoretical  and  experimental 


investigations.  Primary  attention  is  paid  to  the  fundamental  Blausius- 
Mangler  boundary  layer  solution  for  the  sharp  cone,  the  "weak  interaction" 

formulation  of  Probstein  including  transverse  curvature    ,  and  the 

3 
viscous  layer  analysis  of  Cheng  .   The  blunt  cone  results  are  presented 

in  terms  of  parameters  utilized  in  the  local  similarity  analysis  of 

Lees  ,  and  the  thin  shock  layer  formulation  of  Cheng  . 

Becuase  flow  conditions  vary  from  continuum  to  free  molecular, 

theoretical  modeling  for  the  rarefied  conditions  of  the  present 

experiment  is  very  complicated.  The  flow  over  a  significant  portion 

of  a  sharp  body  in  a  rarefied  gas  is  characterized  by  the  "merging"  of 

the  shock  and  the  boundary  layer  near  the  surface.  One  of  the  more 

significant  characteristics  of  the  so-called  "merged  region",  which 

marks  the  transition  from  strong  or  weak  interaction  to  free  molecular 

6 
flow,  is  the  sharply  defined  thickening  of  the  shock  wave.  Feik,  et.al. 

have  investigated  this  phenomenon  on  sharp  cones  of  varying  half  angle 

under  conditions  identical  to  those  of  the  present  study.  Using  a 

radiating  hot  wire  to  probe  the  flow-field  about  the  model,  Feik 

determined  that  for  the  5°  cone  merging  begins  at  a  value  of  V   of 

00 

0.15  to  0.18,  for  the  10°  cone  at  a  value  of  V   between  0.17  and 

00 

0.20,  and  for  the  15°  cone  at  a  V   of  0.20  and  0.22.  At  the  most 

00 

rarefied  conditions   investigated  in  the  present  experiment,   approximately 
2.3"  or  27$  of  the  5°    cone  was  within  the  merged  region,  while  1.0"  or 
33$>  of  the  15°    cone  experienced  merged  flow.     Although  no  3°    cone  data 
is  available,   the  merged  flow  region  is  most  likely  of  considerable 
extent.     As  the   cone   half  angle   is   increased,    it   is  expected  that 
the  boundary  layer  displacement  effect  becomes  more  and  more  dominant 
in  shortening  the  merged  region.     However,   extrapolation  of  Feik's 
data  would  lead  one  to  the  conclusion  that  a  significant  portion  of 
the  20°    and  30°    cone  models  were  immersed  in  merged  flow  during  the 
present  investigation. 


Experimental  measurements  of  heat  transfer  at  hypersonic   speeds 
and  Reynolds  numbers  low  enough  to  produce   "merged  flow"  are  rather 
limited.     Most  of  these   studies  have  been  conducted  in  shock  tunnels 
or  hot  shot  facilities  where  run  times  are  extremely  short.      Shock 
tunnel  heat  transfer  measurements  have  been  performed  on  sharp  and 
blunt  cones  by  Wilkinson  and  Harrington     who  tested  cones  of  6^     and 
9^    at  Mach  numbers  from  13-5  to  22  and  Reynolds  numbers   from  4,000  to 

19,000  per  inch.      Their  sharp  cone  results  followed  the   viscous  layer 

3 
prediction  of  Cheng     at  the   lower  Reynolds  numbers  (w  <  50)  and  tended 

2 
toward  the   transverse  curvature  predictions  of  Probstein     at  the 

higher  Reynolds  numbers.      Their  blunt  cone  results   showed  heat  transfer 

as  much  as  75$  greater  than  had  been  predicted  by  Lees    .     Shock  tunnel 

results  at  lower  Reynolds  numbers  ( .5  <  W  <  10)  were  obtained  by  Waldren 

for  5°  ,   10°    and  20°    sharp  cones  with  Mach  numbers  from  19-24  and  Reynolds 

numbers  from  150-4,000  per  inch.     Waldren 's  results  agreed  well  with 

the  viscous  layer  analysis  of  Cheng.     Other  shock  tunnel  data  for 

9 
sharp  9°    cones  has  been  obtained  by  Burke  and  Dowling     and  for  5°    sharp 

10 

and  blunt  cones  by  Witliff  and  Wilson  .  When  the  various  shock 

tunnel  results  are  compared,  they  exhibit  significant  scatter  which 


is  as  large  as  5O70  even  for  the  higher  Reynolds  numbers.  Griffith 
and  Lewis  '  have  obtained  hot  shot  data  for  5°  spherically  blunted 
cones  at  Mach  numbers  near  20  and  Reynolds  numbers  between  8,000  and 
15,000  per  inch.  Their  results  showed  generally  good  agreement 

with  Lees'  laminar  theory,  and  with  a  slight  modification  of  the 

5 
parameters  suggested  by  Cheng  ,  provided  excellent  correlation  of 

the  data. 

A  method  of  avoiding  the    short  run  times  and  the  associated 

thin-film  thermometry,  which  was  utilized  in  the  investigations  cited 

above,   is  the  well-known  thin -wall  technique.     This  method  has  been 

applied  successfully  at  the  Gas  Dynamics  Laboratory  for  the  measurement 


of  heat  transfer  rates  in  its  continous  facilities.  To  date  there 
have  been  very  few  heat  transfer  investigations  under  hypersonic 
rarefied  flow  conditions  in  continuous  facilities.  Previous  measure- 
ments under  conditions  similar  to  the  present  study  have  been  limited 

12     0  13 

to  those  performed  by  Feik   on  y    sharp  and  blunt  cones,  by  Vas  "  on 

Ik 

10°  sharp  and  blunt  cones,  and  by  Gregorek   on  10°  and  25°  blunt  cones. 

The  investigation  reported  herein  therefore  provides  one  of  the  first 
opportunities  for  comparison  of  data  from  a  continuous  facility  with 
that  taken  in  the  short  duration  facilities.   It  also  provides  the  only 
experimental  heat  transfer  data  available  for  3°  and  30°  cones  under 
the  present  hypersonic  flow  conditions. 

In  the  following  section  a  brief  review  of  the  more  pertinent 
theoretical  studies  is  presented.  In  section  IV  comparisons  are  made 
between  the  theoretical  predictions  and  the  results  which  were  actually 
obtained . 


II 

REVIEW  OF  THEORETICAL  STUDIES 

This  section  contains  selected  theoretical  treatments  of  heat 
transfer  to  pointed  and  blunt  cones  in  viscous  hypersonic  flow. 
Although  no  attempt  is  made  to  survey  the  extensive  literature  on 
the  subject,  it  is  hoped  that  the  analyses  considered  are  representative 
of  the  fundamental  and  more  easily  applied  theories  of  the  classical 
boundary  layer,  weak  interaction,  and  viscous  flow  regimes.  Comparison 
of  the  present  experimental  results  to  the  theories  presented  below 
is  postponed  until  Section  IV. 

A.    In vise  id  Cone  Solution 

Although  the  present  study  concerns  the  viscous  phenomenon  of 
heat  transfer  to  cones  through  a  laminar  boundary  layer,  a  knowledge 
of  the  gas  properties  of  the  inviscid  flow  region  between  the  shock 
and  boundary  layer  is  of  considerable  importance.  In  several 
theoretical  studies  described  below,  the  heat  transfer  at  the  body 
surface  is  predicted  in  terms  of  the  theoretical  inviscid  flow  parameters. 
The  flow  field  in  steady,  supersonic  flow  of  an  ideal,  compressible  fluid 

about  a  circular  cone  with  an  attached  shock  wave  was  treated  by  Taylor 

15 
and  Maccoll  '  in  1933*  If  "the  flow  is  isentropic  and  irrotational, 

the  second  order,  non -linear  differential  equation  which  describes 

the  flows  can  be  expressed  in  terms  of  a  single  velocity  component 

from  which  all  other  variables  may  be  obtained.  The  so-called  Taylor- 

Maccoll  equation  can  be  easily  integrated  to  yield  a  solution  in 

hypersonic  flow  if  one  also  adopts  the  small  density  assumption, 

P  /p  «  1  . 

00 

Inviscid  cone  values  considered  in  the  present  study  were  taken 
from  tabulations  by  Sims  (personal  communication)  who  had  performed 
the  integration  in  the  necessary  Mach  number  range.   Inviscid  temperature, 
pressure,  and  density  ratios  were  assumed  functions  of  M  and  Q       alone, 

co         c 

while  the  theoretical  values  of  Reynolds  number  at  the  cone  surface  was 
assumed  a  function  of  M   .  9       .    P   .  and  T 


B.    Blausius -Mangier  Boundary  Layer  Solution 

The  laminar  boundary  layer  on  a  cone  has  been  treated  extensively 
(For  a  review  of  theoretical  results  see  Ref.  l6,   YJ ,   and  18).   The 
major  difficulty  in  solving  the  governing  differential  equation  arises 
from  the  fact  that  the  fluid  properties  (u,  k,  p,  c  )  vary  with  pressure 

Jtr 

and  temperature.     As  a  consequence  certain  approximations  concerning 
the  viscosity-temperature  law  in  the  boundary  layer  must  be  made   or 
an  alternate  empirical  approach  must  be  taken.     A  well-known  approximate 
solution  with  general  applicability  has  been  derived  by  means  of  a 
Crocco  transformation  of  the   compressible  laminar  boundary  layer 
equations.      Through  a  change   of  independent  variables,   the  equations 
of  continuity,  motion  and  energy  are  reduced  to  two  equations  for  the 
unknown  enthalpy  and  shear  stress.     For  the  case  where  the  viscosity 
law  (u  ~  Tw)  may  be  considered  linear  (w  =  1  or  pu  =  constant),  the 
differential  equation  may  be  transformed  to  the  Blausius  equation 
whose  solution  yields  the  result: 

C.     yrRe"      =       0.661+12  II-l 

f       v       x 

It  may  also  be  shown  that 

Nu  =  0.332      \/Re~      Vpt~  II-2 

Since    St  =  Nu/RePr 
one  may  write: 

£i'-2/3 
Vie 


O.332  pr 
St  =   DD         11-3 


This  solution  for  the  nondimensional  heat  transfer  to  a  flat  plate 


with  zero  pressure  gradient   in  the  boundary  layer  can  be  applied  to 
the  cone  where  the  Mangier  transformation   is  valid.     With  the   restriction 
that  the  boundary  layer  thickness  remain  small  compared  to     r  (x),   Mangier 
has   shown  that  a  direct  transformation  of  the  boundary-layer  equations 
for  axisymmetric  flow  to  those  of  two-dimensional  flow  is  possible  for 
a  perfect  gas  with  constant  specific  heats.     The  heat  transfer  at 
certain  distances  along  the  slant  surface  from  the  tip  is  shown  to  be 
y3  times  the  heat  transfer  at  the   corresponding  flat  plate  distances 
from  the  leading  edge . 

Thus  the  Blausius -Mangier  solution  for  heat  transfer  on  a  cone 
may  be  expressed: 


St 


c 


V3~(0.332)Pr"2/3  llJ+ 


(Re     J 
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It  has  been  found  that  a  better  estimate  of  compressible  heat 

transfer  may  be  obtained  by  a  modification  of  the  above  equation 

18 
according  to  the  reference  enthalpy  method  of  Eckert   .  The  basic 

assumption  of  this  empirical  method  is  that  heat  transfer  can  be 

calculated  with  incompressible  flow  relations  provided  all  temperature 

dependent  gas  properties  are  evaluated  at  an  appropriate  reference 

enthalpy.  According  to  Eckert,  the  Blausius -Mangier  equation  should 

be  altered  so  that : 


st  .  yr(o.332)Pr-2/3    /V. 


VRe 


x.c  e  e 


where  the  starred  quantities  are  to  be  calculated  at  the  reference 

enthalpy.  However,  instead  of  using  Eckert 's  reference  enthalpy  or 

19 
the  wall  temperature  as  a  reference  as  suggested  by  Chapman  and  Rubesin  , 

one  may  proceed  in  an  alternative  manner.  Under  the  assumptions  that 


the  pressure  is  constant  across  the  boundary  layer  and  that  the  "body 
is  slender,  one  may  write: 


/?*"*     /^Te       /^ 
V  Pe^e  =  V  p^e   S  V  ^T* 

The  choice  of  a  reference  temperature  for  this  temperature  - 

viscosity  relationship  may  he  determined  according  to  the  expression 

5 

proposed  by  Cheng   : 

T*     =     T       [1  +  3(Tw/T0)]/6  II-T 


For  a  locally  similar  boundary  layer  flow  around  a  flat  plate  Cheng 
applied  the  Crocco  integral  to  show  that     T       is  the  average  temperature 

TV 

in  the  boundary  layer.  Thus  one  has  the  modified  Blausius -Mangier 
equation  which  is  applied  in  the  present  study: 


-2/3   r* 


Stc  =  VT" (0.332)Pr"  n    J  — II-8 


where 


C 


*      U  T 
oo  * 

Although  a  solution  of  this  form  is  important  because  of  its  wide 
applicability  in  the  continuum  range >   it  also  constitutes  the  zeroth- 
order  solution  for  regimes  of  hypersonic  flow  where  conditions  vary 
considerably  from  those  of  the  thin  boundary  layer.  The  present  work 
verifies  that  by  the  choice  of  an  appropriate  reference  temperature  the 
validity  of  this  solution  may  extend  into  gas  dynamics  regimes  where 
the  assumptions  of  its  derivation  are  inapplicable. 


C.    Boundary  Layer  Displacement  Effect  and  Transverse  Curvature 

l)  Weak  Interaction 

The  laminar  hypersonic  boundary  layer  on  a  cone  is  characterized 
by  a  relatively  thick  curved  boundary  layer  which  produces  outward 
deflections  on  the  stream  surfaces.  This  effect  can  be  large  enough 
to  significantly  change  the  "effective"  shape  of  the  body  and  perturb 
the  inviscid  conical  flow  field.   The  perturbations,  which  result 
from  "induced"  pressure  changes,  are  transmitted  back  into  the  viscous 
boundary  layer  and  affect  its  characteristics.  The  phenomenon  was 
first  recognized  by  Probstein  who  authored  the  definitive  work  in  the 
the  weak  interaction  regime.  Here  the  influence  of  the  shock  wave 
on  the  almost  uniform  flow  results  in  perturbations  that  are  small 
compared  to  those  upstream  where  the  shock  is  highly  curved. 

Probstein' s  model  of  the  flow  consists  of  a  distinct  viscous 
layer  of  hot  gases  near  the  surface  separated  by  a  much  cooler  zone  of 
rotational  inviscid  flow.  The  analysis  is  governed  by  the  Prandtl 
boundary  layer  equation  neglecting  the  influence  of  the  immediate  nose 
region  and  slip  effects.  The  induced  pressure  field  is  considered  to 
be  specified  by  the  induced  flow  inclination  angle  at  the  edge  of  the 

viscous  layer.  Representation  of  this  pressure  field  is  obtained  (in 

do* 
a  Taylor  series)  in  powers  of  the  total  flow  deflection  angle,  - —  . 

The  coefficients  of  the  expansion,  which  are  in  the  form 

1  rd    P/fto,  1  *2    P/Poo 

7 —     [— : 3  and     7 —     [ r-1  II -9 

vj^       d  e     e^c  V£        a  e2    e=ec 

may  be  expressed  as  functions  of    7,  M       and  the  hypersonic   similarity 

00 

parameters  K  and  K  •  It  is  assumed  that  the  tangent  cone  approximation 

and  hypersonic  similarity  apply  so  that  K  must  be  _>  1  and  0     <  20°  . 

do*  ° 

Under  these  conditions  the  value  of  - —  induced  by  the  viscous  flow 

dx 
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is  shown  to  depend  on  the  pressure  distribution  p(x,S).  Upon  adoption 

of  a  linear  viscosity  relation  after  the  form  of  Chapman -Rube sin,  the 

induced  pressure  field  is  expressed  in  terms  of  the  interaction 

parameter  X   .  The  velocity  and  temperature  profiles  of  the  inter- 
c 

action  viscous  field  may  then  be  expressed  asymptotically  in  ascending 

powers  of  X   .  Numerical  solutions  for  the  first-order  corrections 
c 

to  the  velocity  and  temperature  distribution  for  Pr  =  1  were 
obtained  by  substitution  of  the  asymptotic  expressions  into  the 
boundary  layer  equation  with  the  outer  boundary  condition  determined 
by  a  term  for  term  matching  with  the  expansion  for  the  external  viscous 
flow.  The  resultant  first -order  correction  for  the  boundary  layer 
displacement  effect  is  expressed 


*'\ 


=        [-0.350    +    0.111    C^/rp     )     +    0.429   7     " 


\  _  11-10 

0.028  M     ][(d     F  (k)  i|r    )/7M     ] 

The  zero  order  solution  is  represented  by     a       which  is  the  classical 
Blausius -Mangier  solution  discussed  above. 

2)     The  Transverse  Curvature  Effect 

The  laminar  boundary  layer  of  an  axisymmetric  body  differs  from 
that  of  a  two-dimensional  body  because  it  must  spread  circumferentially 

as  it  grows  with  distance  along  the  surface.  An  investigation  of  this 

2 
effect  by  Probstein  and  Elliot  began  with  the  premise  that  the  rate 

at  which  the  body  circumference  changes  with  length  is  the  axisymmetric 
geometrical  factor  which  determines  the  characteristics  of  the  viscous 
layer c  It  was  found  that  the  transverse  curvature  effect  is  characterized 
by  the  parameter  /\/r   (x),  which  is  essentially  the  ratio  of  boundary- 
layer  thickness  to  body  radius.  For  an  insulated  cone  of  half -angle  9 
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in  supersonic  flow: 

0.99^  +  0.275  M  2        r~n 

t^oW)     =   SI -"VST*  n-n 

M  c  x 

This  expression  is   intended  as  a  criterion  for  determining  the 
extent  of  the  nose   region  where     A/r  (x)  »  1     and  the  downstream 
region  where     A  A  (x)      is  on  the  order  of  or  less  than  unity. 

Probstein  and  Elliot  treat  the  transverse  curvature  effect  in 

the  region  where     a  A"  (x)  <  1   .     First,  the  Busemann  and  Crocco  integrals 

of  the  two-dimensional  energy  equation  for     Pr  =  1     are   shown  to  be 

valid  for  axisymmetric  flow  in  which  transverse  curvature  effects  are 

considered.     Then  the  "boundary-layer  equations  are  reduced  to  an  almost 

two-dimensional  form  through  a  modified  Mangier  transformation,  and  a 

solution  is  obtained  for  the   region  where     a A  (x)  <  1   .     For  the  cone 

o 

with  zero  pressure  gradient  in  a  perfect  gas  having  constant  specific 
heats  and  a  Prandtl  number  of  unity,  the  heat  flux  may  be  written: 

=  [0.517  +  0.913  (  w/t  )  +  0.121  (7-I )M     ] 

Si  °  c  11-12 

[VT tan  ac^T] 
x  c 

This  equation  is  also  expressed  in  terms  of  the  classical  boundary- 
layer  results,   q^  ,  which  one  obtains  using  the  Mangier  formulation. 
It  is  important  to  note  that  the  corrections  obtained  for  the 

boundary-layer  displacement  effect  were  also   in  inverse  powers  of 

l/2 
(Re        )        .     Probstein  recognized  that  the  principle  of  linear  super- 

position  applies  and  that  the  transverse  curvature  and  induced  pressure 

effects  are  directly  additive.     By  combining  the  corresponding  corrections, 

Equations  11-10  and  11-12,  the  total  heat  transfer     may  be  expressed  as 
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Stc  =   -/3     (0.332)  Pr"2/3 1    ^-)      {1  +   [-0.350  +  0.111  (~J 

X   c  c 

2  —  2 

+  (0.429  7    -  0.028)M     ][d     F,(k)  \Ir  /7M     ]      + 

c         c     1  c       c 


+   [0.517  +  0.913(Tw/tJ   +   0.121(7-1)M     ] 


11-13 


1  .     _C_v 

L/3"tane  .  Re    '    H 

*-  c 

The  transverse  curvature   correction   is   larger  than  that  for  induced 
pressure  over  the  entire  cone  half -angle  range.     For  the   present 
test  conditions  and  the  3°    cone  the  transverse  curvature  correction 
is  about  50$  larger  while  for  the  20°    cone  it   is  about  20$  larger. 

D.  Solomon-Loeb     Pohlbausen  Solution  for  Hypersonic  Boundary  Layers 

Since  the  publication  of  Probstein's  papers  on  the   induced 

pressure  and  transverse  curvature  effects,   several  approaches  to  the 

20 

same  or  similar  problems  have  been  treated  by  the  Pohlhausen  method 

Application  of  the  Pohlhausen  technique  involves  an  integration  of 
the  boundary-layer  momentum  equation  to  yield  an  ordinary  differential 
equation  for  the  variation  of  one  free  parameter  along  the  body  surface. 
The  solution  of  this  equation  completely  specifies  the  approximate 
velocity  profile  at  each  point  along  the  body  surface,  and  all  flow 
properties  of  interest  may  be  determined. 

A  recent  Pohlhausen  solution  for  conical  flows  has  been  obtained 

21 
by  Solomon  and  Loeb   who  consider  the  same  problem  as  Probstein's: 

the  approximate  analysis  of  the  self -induced  pressure  on  slender 

bodies  in  the  presence  of  transverse  curvature  effects.  The  only 

major  difference  in  the  formulation  of  the  problem  is  that  Probstein 

assumes  hypersonic  similarity  and  procedes  by  the  series  expansion 

method,  while  the  authors  of  the  more  recent  study  utilize  the 
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Pohlhausen  technique.     Restrictions  on  the  analysis  and  the  other 
assumptions  are  very  similar  (See  Sec.   C).      The   flow  field   is  assumed 
to  consist  of  two  parts:     A  viscous  boundary  layer  adjacent  to   the 
body  surface  and  a  purely  inviscid  flow  between  the  boundary  layer 
and  the   shock.      In  the  boundary  layer     Pr     and     C       are   considered 

XT 

constant  and  a  linear  viscosity-temperature  law  is  assumed.   The 

inviscid  Mach  number,  M  ,  is  considered  large  enough  so  that  H  /h  »  1. 

'  e'  e  e 

The  latter  assumption  allows  the  authors  to  adopt  the  hypersonic 

2 
approximation  which  essentially  specifies  that  u  =  2H 

e     e 

The  governing  momentum  integral  equation  obtained  by  the  Pohlhausen 

technique  is  convenient  because  the  inviscid  flow  enters  only  through 

a  variation  of  P  .  The  resultant  pressure  dependence  suggests  the 
e 

adoption  of  an  empirical  expression  for  the  inviscid  pressure  field 
in  order  that  the  induced  pressure  problem  may  be  studied  in  conjunction 
with  the  transverse  curvature  effect.  This  induced  pressure  problem 
can  only  be  solved  where  the  boundary  layer  and  the  inviscid  flow  are 
solved  simultaneously.  Thus,  the  Pohlhausen  solution  of  the  boundary 
layer  is  constructed  to  allow  for  a  significant  alteration  in  the 
inviscid  flow  field  due  to  the  displacement  of  the  boundary  layer. 
To  solve  for  the  inviscid  flow  field,  an  "effective  body"  is  assumed 
to  exist  with  the  same  inviscid  field  as  the  original  body  with  all 
viscous  effects  neglected.  The  surface  pressure  on  the  "effective 
body"  is  described  by  a  tangent  cone  formula  in  the  following  form: 

p  2 

e/p   =  a)  K  +u),     where    K  =  M  tan  0  II-14 


:/p    =  0)  K  +  U)n 
'  oo        O         1 


00 


In  this  manner  the  local  inviscid  pressure  on  the  "surface  of  the 
effective  body"  is  related  to  the  undisturbed  flow  ahead  of  the  body 
through  conical  flow  theory,  but  using  the  local  surface  inclination 
as  the  effective  cone  angle. 
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Solutions  of  the  governing  momentum  integral  equation  with  the 

inviscid  pressure   so  defined  are  obtained  by  Solomon-Loeb  for  the 

limiting  cases  of  very  small  and  very  large  transverse  curvature  effects. 

The  limiting  cases  are  defined  in  terms  of  a  variable  which  contains 

all  of  the  relevant  parameters: 

k 
0)     tan  G  Re 
—  o  C      X,oo 

X     =     ~ £*-  11-15 

(7-1)  C     M 

*         00 

Since  the  transverse  curvature  effect   increases  with  decreasing  Reynolds 

numbers,    increasing  Mach  number,  and  decreasing  semi-vertex  angle,    it 

is  quite  natural  that  the  heat  transfer  solution  for  large  transverse 

curvature  be  expressed  as  a  limiting  solution  as     X  -  0   .     Neglecting 

terms  of     0([log  log  =]/log  X)     the  heat  transfer  may  be  expressed: 

A. 


2(u     A  _  H         HP   -R 

V     7-1     '  c  HePr  ^   H^   ; 


C_         =     : as     X     -    0  II-16 

x  log  f 

For  the  opposite  limit  of  small  transverse  curvature  effects,  the 
heat  transfer  solution  is  given  in  decreasing  negative  powers  of 
X     and  the  classical  boundary  layer  solution  is  approached  as     X  -» °° 
For  this  case  the  equation  may  be  written 


H  Pr-H      co     tan     9      JUto, 
,    e         w.     o c    v 1       1 


2/     1    1 


+    h  (Tg-)  t3  (^A +  n-17 


f  (1  ♦  3£^0)]  4-  )     as     X  — 
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E.  The   Viscous -Layer  Regime  as  Analysed  by  Cheng 

Due  to  the  design  alterations   in  the  20°  and  30°    cones,    it  was 
possible  to  obtain  heat  transfer  data  at   stations   in  the  vicinity  of 
the  leading  edge  (0.1  inch  <  x)   of  these  models.     For  the   lower  Reynolds 
number  conditions  of  the  present  study,   the  flow  about  the  forward 
regions  of  the  models  denoted  above  may  be   classified  as  belonging  to 
the   "viscous -layer  regime",   because  here  the  viscous  layer  becomes 
the  dominant  fraction  of  the  shock  layer.     Departures  from  the  predictions 
of  thin  boundary-layer  concept  in  this  flow  regime  become  more  and  more 
evident  as  rarefaction  is  increased  (See  Ref.   8).     To  the  author's 

knowledge,   the  only  solution  applicable  to  flows  over  pointed  cones 

3 
where  the  shock  layer  is  fully  viscous  has  been  formulated  by  Cheng   . 

The  outstanding  features  of  the  portions  of  Cheng's  analysis  which 

are  applicable  to  the  present    investigation  are  discussed  briefly 

below. 

According  to  Cheng's  flow  model  for  the  viscous  layer  regime, 

the  flow  field  around  a  body  in  rarefied  flow  can  be  divided  conveniently 

into  two  regions:     the  shock  layer  and  the  shock-transition  zones. 

The  compression  ratio  across  the  shock-transition  zone   is  assumed  to 

be  high,  and  both  regions  are  assumed  to  be  thin  in  comparison  to  a 

typical  body  dimension.     The  analysis  begins  with  a  consideration  of 

the  Navier-Stokes  equations  governing  compressible,   heat  conducting 

flow.     By  neglecting  transverse  curvature  and  retaining  only  the 

highest  order  normal  derivatives  in  the  transport  terms,  Cheng  obtains 

a  parabolic   system  of  equations   similar  to   the  boundary -layer  equations 

of  compressible  flow.     In  the  formulation  of  the  problem  for  the  thin 

shock  layer  the  non-slip  boundary  condition  is  adopted  at  the  body 

surface.     This  boundary  condition  is  justified  by  a  manipulation  of 

the  momentum  and  energy  equations  of  the  Navier-Stokes  equations  to 

show  that  the  phenomenon  of  velocity  slip  and  temperature  jump  can 
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be  neglected.  Other  assumptions  of  the  analysis  involve  the  neglect 
of  transverse  curvature,  the  boundary  layer  displacement  effect  and 
shock  curvature  in  comparison  to  the  more  important  slip-like  velocity 
and  enthalpy  changes  across  the  shock.  The  validity  of  the  analysis 
is  limited  to  thin  shock  layers  over  highly  cooled  surfaces  and  non- 
slender  cones  where  5/r  «  1  .  The  gas  is  assumed  to  be  perfect 
with  constant  specific  heats  and  unit  Prandtl  number.  Under  this 
framework  and  with  the  adoption  of  a  linear  viscosity  law,  the 
analysis  is  reduced  to  the  solution  of  a  single  parabolic  differential 
equation  governing  the  velocity: 

3(u%)    _3(u%)       32(u%) 

w  ~w~  - 2  *  ~W~  +  k<  ~tf-  IJ-18 

with  the  boundary  conditions  (f  =  1  at  shock): 
u(0)   =  0 


u(l)-l+f 


.,-2       -,  11-19 

o(u  /  ) 


d\|r 


J     *  =1 


Since  the  equation  is  of  the  parabolic  type,  forward  integration  is 
possible  and  a  nonsimilar  solution  may  be  obtained  by  the  method  of 
finite  differences.  Comparison  of  the  resultant  solution  to  the 
Blausius  profile  reveals  close  agreement  down  to  a  value  of  W  ~  3  > 
even  though  the  shock  layer  is  fully  viscous  at  that  point.  The 
comparison  reveals  that  departure  from  the  predicted  values  of  heat 
transfer  by  the  thin  boundary  layer  theory  should  first  occur  at 
W  =   3  and  become  significant  at  W  =  1  . 

It  is  important  to  note  that  the  effects  of  gas  rarefaction  and 
surface  incidence  are  all  contained  in,  and  can  therefore  be  correlated 
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by,  the  variable  W  . 

e   p     u     x     u     T 

00        oo  ,      0         *  » 

W     =     ( — )   sec  e 

LL  >        I  S 

0*0 
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2      y  C     cos  0 
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Cheng  presents  his  heat  transfer  prediction  in  terms  of  the  parameters: 

CHoo/  vs       W  11-21 

H>°°/sin  e 
s 

The  numerical  solution  presented  is  valid  in  the   region     0  <    W  <     10   , 

and  it  approaches  the   classical  free -molecular  and  boundary  layer  results 

at  the  extremes  of  the  region.     The  solution  first  begins  to  depart 

from  the  classical  boundary  layer  theory  at     W  ~  3     an(i     C  approaches 

H,°° 

the  value  of  sin  6       as  W  -♦  0  as  required  by  free -molecular -flow 
theory  for  unit  accommodation  coefficients. 

Cheng's  formulation  is  especially  valuable  in  the  analysis  of  the 
present  data  because  ail  of  the  experimental  variables  are  contained 
in  one  convenient  correlation  parameter. 

F.    Lees'  Local  Similarity  Analysis 

The  concept  of  assuming  a  local  similarity  in  the  boundary  layer 

as  a  means  of  estimating  heat  transfer  to  blunt  bodies  at  high  speeds 

k 
was  first  proposed  by  Lester  Lees  .  By  using  the  Levy  transformation 

to  transform  the  momentum  and  energy  equations  into  a  more  convenient 

"low  speed"  form,  he  obtained  locally  similar  solutions  to  the 

boundary  layer  equations.     Lees  observed  that  the  enthalpy  distribution 

in  the  boundary  layer  is  rather  insensitive  to  pressure  gradient .      In 

effect  he  found  that  the  nondimensionalized  enthalpy  gradient  at  the 

surface  is  represented  very  accurately  by  the   classical  zero-pressure 

gradient  value,  and  that  the   surface  heat  transfer  rate  can  be  obtained 
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in  terms  of  the  pressure   distribution. 

Applications  of  the   result  were  obtained  for  the   special  cases 
of  the  hemisphere  and  the  right  circular  cone   capped  by  a  spherical 
segment.     The  modified  Newtonian  pressure  distribution  in  the  form: 

2  2  2 

P/n       =     cos  e      (1  +  7     M     )   sin  G  11-22 

■Cq  00  00 

was  adopted  as  a  means  of  calculating  the  heat  transfer  distributions. 

For  the  hemisphere  and  the  leading  edge  of  the   spherically  blunted 

cone  Lees'   result  may  be  written: 

2  1  2  2 

<1  (1/2)9   sin  e    [cos  Q  +(   jy  M      )sin  0] 

W  oo    co  „ 

11-23 


(q)  r°°  2r        2  ,1/  2.  2   - 

K^w'  J      sin  9   [cos  9  +  (    h  M     )   sin  0]d0 

O  O  oo    oo     ' 

The  modified  Newtonian  pressure  distribution  is  also  assumed  for  the 
conical  skirt  where  all  physical  quantities  are  regarded  as  constant. 
The  result  for  the  heat  transfer  distribution  for  x/r  ^  cot  6       may 
be  written 

%  X/ro 

*  A(«  )  ,  w,  n-2^ 


{^l  c     WO  +  (xAo)3l  1/2 


A(0  )  and  B(S  )  are  functions  of  y ,  M  and  0  (Ref.  4).  At  the 

C  C  oo  C 

junction  of  the  nose  and  conical  skirt  the  two  expressions  for     q/(a    ) 
are  identical.     Lees'   result  predicts  for  blunted  cones  a  heat  transfer 
rate  that  is  smaller  than  that  of  an  equivalent  sharp  cone  for  several 
nose  radii  aft  of  the  nose   -  skirt  junction. 


G.  Combined  Leading  Edge  Bluntneos  and  Boundary  Layer  Displacement   - 

Cheng's  Thin  Shock  Layer  Theory 

Within  the  formalism  of  the  classical  boundary  layer  theory  and 

3 
the  hypersonic  small  disturbance  theory,  Cheng     has  theoretically 
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treated  hypersonic  flows  over  thin  bodies  subject  to  the  combined 
effect  of  leading  edge  bluntness  and  boundary  layer  displacement  for 
two  dimensional  and  axisymmetric  flows.  Although  separate  treatments 
of  either  boundary-layer  displacement  or  bluntness  alone  received 
considerable  attention  during  the  1950* s>   Cheng  was  the  first  to 
attack  the  combined  problem  satisfactorily.  His  approach  is  charac- 
terized by  two  basic  concepts:  a  local  similarity  in  the  boundary 
layer  and  a  detached  shock  in  the  outer  inviscid  flow.  His  basic 
flow  model  consists  of  a  low  density  entropy  layer  situated  between 
the  laminar  boundary  layer  and  the  shock  layer.  This  inviscid 
region  termed  the  "entropy  layer"  is  so  designated  because  it 
consists  of  gas  particles  which  have  come  through  the  stronger,  forward 
portion  of  the  shock.  The  terminology  is  appropriate  since  the  specific 
entropy  of  these  particles  is  much  higher  than  that  of  particles 
passing  through  the  comparatively  straight  shock  downstream  of  the 
nose  region  and  farther  away  from  the  body. 

The  domain  of  validity  of  Cheng's  formulation  is  limited  by 
several  restrictions  which  govern  the  accuracy  of  the  analysis  in 
the  invisc id-flow  regime  and  the  boundary  layer.  Underlying  the 
analysis  are  the  requirements  that  l/(M  6      )  and  e  =  J -1/7+1  be 
sufficiently  small  in  comparison  to  unity  to  permit  treatment  by  the 

shock  layer  concept.   (For  Mach  numbers  and  flow  deflection  angles 

/  2   2 
such  that  1/M  9       «  1  a  strong  shock  exists).  Through  a  trans- 
formation of  the  boundary  layer  equations,  the  assumption  that  e 
be  small  is  shown  to  be  essential  in  order  to  neglect  the  normal 
pressure  gradient  across  the  boundary  layer.  The  Cheng  formulation 
is  ideally  applied  in  the  limit  of  7—1  because  the  direct  effect 
of  pressure  gradient  is  at  most  of  order  €  .  However,  it  is  argued 
that  the  formulation  is  valid  for  small  e   since  the  leading  approxi- 
mation of  the  solution  to  the  transformed  boundary  layer  equations 
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is  independent  of     e     and  is  governed  by  the  Blausius  equation. 

Cheng's  analysis   is   further  restricted  by  the   simplifying 

assumption  that     u    «  u       at  the  outer  edge  of  the   boundary  layer  and 
e         co 

within  the  entropy  layer  as  well.     He  shows  that  the  error  in  this 
assumption  may  be  expressed: 


u  -u 
,   e     co 

'  u 


<  i1/2)e  2(7_1)/r  n-25 


Cheng  insists  that  this  quantity  be  negligible  in  comparison  to 
unity,  but  at  the  same  time  he  recognizes  that  9       must  not  be  so 
small  that  the  strong  shock  assumption  is  violated.  Satisfaction  of 
the  former  restriction  on  6       not  only  insures  that  the  perturbation 
velocities  will  be  small,  but  also  that  the  analysis  of  the  entropy 
layer  can  be  based  upon  the  fact  that  the  prevailing  density  is  much 
lower  than  that  of  the  free  stream.  As  an  immediate  consequence  the 
pressure  variation  across  the  entropy  layer  can  be  shown  to  be 
negligible.  In  fact  he  demonstrates  that  the  fields  of  pressure 
and  velocity,  as  well  as  the  shock  locations  can  be  consistently 
determined  under  the  restriction  he  sets  forth. 

On  the  formulation  of  the  differential  equation  for  the  location 

of  the  entropy  layer,  Y™(x)  ,  Cheng  considers  the  inner  edge  of  this 

layer  as  an  equivalent  body  surface  to  account  for  the  displacement 

effect.  The  local  pressure  at  the  outer  edge  of  the  entropy  layer 

and  hence  in  the  complete  region  Y,  <  y  <  Y   can  be  calculated 

b  e 

from  the  Busemann  formula  when  the  local  shape  of  the  outer  edge  of 
the  entropy  layer  is  determined: 


2  , 

P.      =     7     P     M  [Y     Y   ']  11-26 

u  co     co     co  e     e 


Since     u(x)  «  u       everywhere  in  the  entropy  layer,  the  continuity  and 
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momentum  equations  are  integrated  twice  to  yield 


[Y     -  Yj7     E     =  constant  11-27 

e         b  d 


A  combination  of  the  two  equations  yields  the  necessary  ordinary 
differential  equation  for  the  entropy  layer  (assuming     j  -»  l): 

[Y     -  Yj[Y     Y  V   =     constant     =     i/^Vek  t  11-28 

e         bee 

Where  the  constant  on  the  right  hand  side  may  be  obtained  from  an 

integral  of  the  momentum  equation  over  a  cross   section  normal  to  the 

•* 
free -stream  direction.      Since     Y,    =  Y     +5      ,    it   is  clear  that  the 

b         w 

manner  in  which  the  boundary  shock  and  the  entropy  layers  interact  is 
contained  in  the  ordinary  differential  equation  given  above. 

For  the  axisymmetric  case  the  necessary  differential  equation 
for     Y  (x)     may  be  expressed: 


e 


-  £        (ZZ')!    =     2  11-29 


where 


Z     =     2  9     Y     /    yek     t  11-30 

c     e 

t,     =     2  Q       S  /Vek     t  11-31 


% 

k     =     P/      Ti  \  I]C"32 

(l/2)p     u  -(rtt'VM 


00  00 


Although  there  results  a  power  law  relationship  for     Z     in  the 
solution  for  small     £     (near  the  leading  edge)  the  asymptotic  behavior 
of  this  equation  for  large     £     yields  a  solution  which  is  expressed 
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in  terms  of  Bessel  functions  of  the  first  and  second  kind.     As  a  result 
the  boundary  layer  pressure   dependence  on  axial  location,   P(£)>    ^s 
characterized  by  the  exponentially  decaying,   oscillatory  behavior  of 
the  Bessel  functions.     The   implication  of  this  result  in  regard  to  the 
theoretical  heat  transfer  predictions  is  obvious.     To  the  leading 
approximation  for  small     e    ,   Cheng  shows  that  the  heat  transfer  may 
be  expressed  explicitly  as  a  function  of  pressure  alone. 

,       _  P/P 

Cw      M^  «  X  °°      •  11-33 

H>°°  t      Ls    p/p   ds1 

V   o  °°    L 

As  a  consequence  of  this  pressure  dependence  a  correspondingly  oscillatory 
behavior  is  manifested  in  the  predicted  heat  transfer  distribution. 

For  the  blunt  cone  Cheng  suggests  correlating  the  heat  transfer 
dependence  in  terms  of  the  parameters . 

2  ,      .i/k  ,     2     -  9c     S  „, 

M      (ek)   '     CTT      e         X+     vs     -== Tl-3>k 

B>     c         t  y€k     t 

The  nature  of  the  oscillatory  solution  may  be  observed  in  Figs.  22-26. 
A  rapid  changeover  to  oscillatory  behavior  is  predicted  for  the   vicinity 
of 

e2    s 

-       =1  11-35 


ek     t 

Cheng  suggests  that  the  oscillation  is  attributable  to  a  repeated 
reflection  or  skipping  of  the  reattaching  shock  layer  as  it  reaches  the 
vicinity  of  the  surface  of  the  conical  afterbody.  However,  he  concedes 
that  the  phenomenon  of  glancing  re impingement  of  the  shock  layer  may 
be  substantially  altered  by  higher  order  calculation  in  e    . 

In  spite  of  the  fact  that  the  several  essentially  restrictive 
assumptions  of  the  analysis  can  nver  be  absolutely  satisfied,  there 
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appears  to  be  a  fairly  extensive  range  of  hypersonic  flows  where  such 
a  "thin  shock  layer"  theory  should  apply.  The  geometries  and  flow 
conditions  of  the  present  investigation  were  such  that  most  restrictions 
were  well  satisfied  although  the  5°  and  30°  cone  angles  approach  the 
lower  and  upper  limits  required  of  9 
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III 
THE  EXPERIMENTAL  METHOD 

A.  Introduction 

The  experimental  method  used  to  determine  the  local  convective 
heat  transfer  to  the  conical  bodies  was  the  well-known  transient 
"thin-wall"  technique.  Sharp  and  blunt  thin-skinned  conical  shells, 
which  varied  in  half  angle  from  3°  "to  30°  >   were  installed  with  copper 
constantan  thermocouples  and  tested  in  the  hot  nitrogen  facilities. 
From  the  simple  heat  balance  for  an  element  of  model  skin  it  can  be 
shown  that  determination  of  the  temperature  rise  with  time  is  necessary 
and  sufficient  for  computing  the  local  heat  transfer  coefficients  at 
the  model  wall.  The  models,  initially  maintained  at  a  low  isothermal 
wall  temperature  (T  /t  =  O.l^),  were  suddenly  exposed  to  the  free 
stream  when  a  pneumatically  operated  shield  was  removed.  The  temperature 
vs.  time  history  of  the  models  was  recorded  on  an  oscillograph  which 
monitored  the  voltage  level  of  each  thermocouple.  The  temperature 
data  was  differentiated  and  digitized  auto-mechanically  in  order  that 
repetitive  calculations  could  be  accomplished  by  computer. 

B.  Facilities  and  Test  Conditions 

Two  nearly  identical  hypersonic  wind  tunnels  were  used  for  the 
present  investigation.  The  majority  of  the  work  was  completed  in  the 
Princeton  Hypersonic  Tunnel  N-3  while  the  Hypersonic  Tunnel  N-5  was 
utilized  for  a  small  portion  of  the  program.  Both  tunnels  are  blow 
down  facilities  with  maximum  running  times  of  25  minutes.  Since  the 
newer  N-5  tunnel  is  an  almost  exact  duplicate  of  the  N-3  facility, 
separate  descriptions  are  redundant  (Fig.  1).  Only  those  details 
having  a  direct  bearing  on  the  present  investigation  will  be  considered, 

since  a  comprehensive  description  of  the  N-3  facility  is  given  by  Vas 

22 

and  Koppenwallner  . 
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The  working  gas   is  high  purity  nitrogen  (contamination  of  oxygen 

less  than  8  parts  per  million)  which   is  heated  by  a  graphite   resistance 

23 
element       to  stagnation  temperatures   sufficiently  high  to  yield  a  static 

temperature  throughout  the  expansion  which  precludes  the  possibility 

of  condensation.     A  wide  range  of  flow  conditions   is  achieved  by 

alternately  exmploying  one  of  two  conical  nozzles  which  expand  the 

flow  into  the  six  inch  diameter  test   section.      The  higher  range  of 

Reynolds  numbers  from  IT, 000-41, 000 /in.  are   obtained  with  the  6°    half- 

angle  conical  nozzle  at  stagnation  pressures  from  2,000  to  4,000  psia 

and  stagnation  temperatures  of  2,000°  K  and  1,700° K  respectively. 

Corresponding  Mach  numbers  range  from  IT  to  19 •     Lower  Reynolds 

numbers  from  7,000  to  l4,000/in.  at  Mach  numbers  from  23  to  27  can 

be  obtained  when  the  expansion  is  accomplished  by  the   10°    half -angle 

conical  nozzle.     The  stagnation  pressures  required  were  2,000,   3>000, 

and  5,000  psia  at  a  constant  stagnation  temperature  of  2,000° K   .     The 

free  stream  parameters  corresponding  to  the  five  stagnation  conditions 

at  which  each  model  was  tested  is  presented  in  Table  I. 

At  4,000  psia  the  diameter  of  the  usable  inviscid  core  of  the 

wind  tunnels  is  about  4"  at  the   nozzle  exit,  decreasing  to  2"  at  about 

twelve   inches  downstream.     All  models  were  aligned  to  the  tunnel 

centerline  with  their  tips  located  a  nominal  4"  from  the  nozzle  exit. 

This  positioning  insured  that  the  entire  conical  surface  was  well 

within  the  inviscid  core.     The  reasonable  quality  of  the  core  section 

at  the  station  of  the  model  tips  is  shown  in  Figure  2.     Since  the 

conical  nozzles  are  only  partially  "contoured"  by  the  growth  of  the 

tunnel  boundary  layer,  a  slight  expansion  process  continues  into  the 

test  section  resulting  in  a  Mach  number  gradient  along  the   centerline 

of  0.3  per  inch  for  the  higher  Mach  number  range.     Corrections  for  this 

gradient  are   incorporated  into  the  computerized  data  reduction  program. 
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In  the   interpretation  of  heat   transfer  data  an  accurate 
determination  of  the  thermodynamic   state   of  the   free   stream  gas   is 
of  critical  importance,   because   it   is  both  customary  and  convenient 
to  normalize  local  heat   flux  by  free   stream  conditions   inferred  from 
the  stagnation  conditions  and  measured  pitot  pressure.      In  addition 
the  nondiraensionalized  heat   flux  data  is  usually  correlated  by 
parameters  which  are   strong  functions  of     Re        ,   the  Reynolds  number 

X  ,00 

based  on  free  stream  conditions.  Consequently  any  error  in  measuring 
the  stagnation  conditions  can  be  compounded  when  results  are  compared 
to  theory.  Some  attention  should  therefore  be  given  to  the  assumptions 
and  measured  quantities  used  to  infer  these  free  stream  conditions  in 
the  hot  nitrogen  facilities. 

Although  dissociation  effects  are  of  no  concern  in  the  N-3 
facility,  the  thermodynamic  properties  of  the  gas  must  be  specified 
with  attention  paid  to  the  vibrational  energy  of  the  diatomic  gas  as 
well  as  the  translational  and  rotational  constituents  of  the  total 
thermal  energy.  Consideration  of  the  former  demands  that  the  freezing 

temperature  for  the  vibrational  mode  of  the  gas  be  determined.  Vas 

22 
and  Koppenwallner   have  shown  that  a  freezing  of  the  vibrational  energy 

occurs  so  close  to  the  throat  that  the  assumption  that  the  flow  is 

vibrationally  frozen  throughout  the  expansion  results  in  only  0.5$ 

error  in  Mach  number.  Acceptance  of  this  approximation  results  in  a 

considerable  simplification,  the  number  of  necessary  parameters  being 

reduced  from  four  to  three:  P  ,  T  ,  and  P^  .  With  knowledge  of  these 

0     0         t 

three  quantities  the  thermodynamic  system  is  completely  specified  under 
the  vibrationally  frozen  assumption  and  all  free  stream  quantities 
of  interest  can  be  calculated. 

Impact  pressures  are  measured  using  a  variable  reluctance 
transducer  calibrated  with  a  silicon  (high  vacuum)  oil  manometer, 
while  the  stagnation  pressure,  P  ,  is  measured  with  a  precision 
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Heise  bourdon  gage.  Since  the  elevated  temperature  of  a  highly- 
pressurized  gas  is  difficult  to  measure  accurately  with  thermocouples, 
the  stagnation  temperature  is  inferred  from  the  mass  flow  rate  of  the 
incoming  gas  through  a  Potter  turbine  flowmeter.  Details  of  this  method 
should  be  scrutinized  since  Stanton  number  and  Reynolds  number  are 
functions  of  stagnation  temperature. 

Inference  of  the  stagnation  temperature  from  the  mass  flow  rate 
is  justified  if  a  number  of  assumptions  are  valid.   If  along  the  gas 
inlet  line  at  any  point  of  interest  (including  the  stagnation  chamber 
and  nozzle  orifice)  one  can  consider  the  flow  to  be  isentropic  and 
sufficiently  uniform  to  warrant  the  one -dimensional  assumption,  the 
same  mass  flow  must  exist  at  any  station.  For  steady  flow: 

m  =  p  uA  =  constant  III-l 

Inserting 

p  =  P/RT   and  u  =  mVtRT 
One  may  write: 


m  =  (P/RT)M   Y7RT  A  III -2 

Utilizing  the  isentropic  relationships  for  total  temperature 
and  pressure: 


T0/T     =      1  +  ^i  M2  III-3 

po/p   -    Id  ♦  ^)7/{7'1) 
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One  obtains; 


P~  A 


*=Wr     ^   {-  2f  111-5 

Vt^        (1  +  ^  M2)(7+1)/(7-D 


Imposing  the  sonic  condition  at  the  throat  and  denoting 

1  , 

r  =  —f — rr-, in-6 


(l+Zf^Vafr-l) 

One  may  write  finally: 


r  P  A„ 


-^r1     Fk  in-7 

The  functional  dependence  of  r  on  pressure  has  been  investigated 
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and  is  presented  by  Vas  and  Harvey  .  Knowing  V;   the  mass  flow  rate,  m  ; 

the  stagnation  pressure,  P  ;  and  the  throat  size,  A  ;  the  stagnation 
temperature  can  be  calculated. 

A  Potter  turbine  flowmeter  was  used  to  measure  the  mass  flow. 
Although  the  turbine  flowmeter  is  ideally  a  volumetric  flowmeter,  it 
becomes  a  mass  flowmeter  upon  specification  of  pressure  and  temperature 
of  the  gas  in  the  meter. 

The  output  of  the  flowmeter  is  related  to  the  mass  flow  as 
follows : 


co  -  m 


n/ap  III -8 


where : 


co  =  frequency  of  turbine  flowmeter  (cps) 

p  =  density  of  gas  in  flowmeter  based  on  T   , 

FM 

a   =  experimentally  determined  flowmeter  constant 
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The  output,  a)  ,    is  displayed  on  a  C.M.C  Universal  Counter 
Timer.  The  density  of  the  gas  in  the  flowmeter  was  determined  from 
data  presented  by  Michels,  et.al.    which  accounted  for  the  non-ideal 
behavior  of  nitrogen  under  high  pressures. 

Calibration  of  the  flowmeter  over  the  pressure  range  from  1,500 
to  5>000  psia  allowed  the  determination  of  the  constant,  a  •   It  was 
found  that  the  Potter  meters  used  were  not  true  volumetric  flowmeters 
but  that  a     decreased  with  increasing  stagnation  pressures  for  each 
of  the  three  meters  used.  Each  meter  showed  a  decrease  of  approximately 
7$  in  a     over  the  pressure  range.  Since  an  error  of  this  magnitude 
in  flowmeter  count  would  have  resulted  in  a  stagnation  temperature 
error  of  about  15$  >  accurate  values  of  a    at  each  stagnation  pressure 
were  determined  to  within  an  accuracy  of  1$,  and  a   =  (X   (P-ovt)  'was 
inserted  into  the  flowmeter  equation: 

tpa* 


U)  = 


—  PI*     77T-T,  m-9 


-/~T  "   v  R  a{Pmh 


From  the  above  equation  working  curves  were  constructed  of 
T  vs.  a)   for  varying  inlet  gas  temperatures,  different  stagnation 
pressures  and  throat  diameters.  Although  the  calibration  mentioned 
above  was  performed  at  room  temperature,  it  was  further  assumed  that 
the  throat  boundary  layer  during  tunnel  operation  is  infinites imally 
small,  since  the  Reynolds  number  based  on  orifice  conditions  is  so 
large . 

The  scatter  in  the  "cold  flow"  calibration  used  to  determine 

effective  throat  diameter  (knowing  P  ,  T  ,  and  oo)  limits  one  to 

an  uncertainty  in  T   of  ±  50° K  at  T  =  2,000° K.   A  more  serious 

o  o 

inaccuracy  in  this  method  of  inferring  stagnation  temperature  arises 
from  the  deposition  in  the  nozzle  orifice  of  carbon  particles 
loosened  from  the  graphite  heater  element.  Over  the  period  of  an 
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individual  tunnel  test  the  effective  throat  diameter  sometimes  changes 
several  ten -thou sanths  of  an  inch  from  the  nominal  0.0320  orifice  employed 
with  the  10°  conical  nozzle.  Since  the  stagnation  temperature  varies 
as  the  square  of  the  throat  diameter,  data  vas  discarded  when  changes  of 
this  order  occurred.  However,  minor  fluctuations  in  throat  error  were 
unavoidable.  The  resultant  variations  could  be  detected  by  a  reference 
thermocouple  on  each  model  which  monitored  heating  rates  at  a  particular 
location  during  each  exposure  to  the  flow.  When  the  tunnels  operate 
at  the  lowest  stagnation  pressure  (2,000  psia),  carbon  tends  to 
deposit  gradually  in  the  orifice  decreasing  its  area.  For  a  constant 
flowmeter  count  there  results  a  continually  decreasing  stagnation 
temperature  even  though  the  power  consumption  of  the  resistance  heater 
remains  constant.  Conversely,  at  higher  stagnation  pressures  the 
previously  deposited  carbon  tends  to  blow  downstream  with  the  corre- 
sponding increase  in  throat  area  causing  an  associated  increase  in 
stagnation  temperature.  The  stagnation  temperature  is  therefore 
uncertain  to  ±  100° K  ,  or  ±  5$  •  An  improvement  in  the  construction 
of  the  graphite  resistance  element  would  considerably  reduce  scatter 
in  future  data  to  be  obtained  in  the  facilities.  Implications 
and  further  observations  concerning  the  problems  mentioned  above 
will  be  considered  in  the  error  analysis,  Section  IH-G. 

Cj The  Experimental  Models 

Typical  blunt  and  sharp  heat  transfer  models  are  shown  in 
Figures  3  and  4.  Each  model  consists  of  a  right -circular  conical 
shell  capped  by  a  spherical  segment.  The  "sharp"  models  were  designed 
with  a  spherical  cap  only  0.004"  in  diameter.  Since  this  dimension 
is  on  the  same  order  as  the  mean  free  path,  one  can  consider  a  mathe- 
matically sharp  cone.  The  "blunt"  cones  are  constructed  with  a 
spherical  segment  whose  diameter  is  0.225",  approximately  100  mean 
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free  paths . 

The  conical  shells  were  constructed  of  electroformed  nickel 
and  instrumented  with  ten  to  twenty  0.005"  copper-constantan  thermo- 
couples which  were  peened  tightly  into  place  at  fairly  uniform  intervals 
along  the  conical  surface.  No  attempt  was  made  to  describe  the  heat 
transfer  to  the  cylindrical  portion  of  the  model,  although  at  least 
two  thermocouples  were  installed  there  so  that  a  tunnel  stall  or 
conduction  from  the  sting  could  be  detected. 

As  a  first  step  in  the  construction  of  a  model,  a  plexiglass 
mandrel  was  smoothly  machined  to  the  required  inner  diameter  of  the 
cone.  Nickel  was  then  electroplated  to  a  thickness  of  about  0.020". 
From  this  rough  surface  a  smooth  conical  shell  was  machined  to  a 
uniform  thickness  of  0.010"  and  to  the  desired  apex  angle.  The  model 
was  then  instrumented  with  the  copper-constantan  thermocouples.  Most 
of  the  thermocouples  were  placed  along  a  particular  ray  of  the  cone. 
However,  four  thermocouples  were  placed  in  the  surface  at  90°  intervals 
of  arc  in  order  that  it  could  be  verified  that  the  angle  of  yaw  in  the 
wind  tunnel  test  section  was  as  close  as  possible  to  zero  degrees.  The 
cylindrical  afterbody  of  the  cone  was  then  attached  to  a  water-cooled 
holder  which  provided  two  leadout  tubes  for  the  thermocouples,  an 
inlet  tube  for  the  nitrogen  cooling  gas,  as  well  as  a  means  for  sting 
mounting  the  model  in  the  wind  tunnel. 

Initially  the  heat  transfer  models  were  designed  with  a  1^"  base 
diameter  and  interchangeable  tips.  Two  of  these  models,  denoted  as 
"Type  I",  are  depicted  in  Figure  3-  Type  I  models  were  designed  with 
interchangeable  tips  for  convenience  and  economy.  A  single  model  of 
this  type  would  suffice  both  as  a  "sharp"  and  a  "blunt"  model  with  a 
single  tunnel  installation,  and  conversion  could  be  accomplished  in 
minutes  without  affecting  the  alignment  of  the  model  to  the  flow. 
The  versatility  of  the  design  is  of  significance  in  that  the  effects 
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of  bluntness  could  not  be  obscured  by  a  variance  in  model  skin 
thickness  or  a  change  in  tunnel  alignment  resulting  from  the  installa- 
tion of  a  separate  sharp  or  blunt  model  of  the  same  apex  angle. 

Upon  the  reduction  of  the  data  from  the  first  Type  I  model  to 
be  tested,  a  disadvantage  of  this  construction  was  revealed  which 
led  to  a  modification  in  the  design  of  several  models.  The  solid 
steel  tip  and  its  mount  inside  the  nickel  shell  constituted  a  sub- 
stantial heat  sink  which  resulted  in  considerable  conduction  errors 
in  the  tip  region.  As  a  result  a  significant  portion  of  the  data 
from  the  more  rarefied  region  of  the  shell  was  discarded. 

The  conduction  to  the  solid  tip  and  the  necessity  of  reducing 
the  model  base  diameter  resulted  in  the  construction  of  the  Type  II 
models.  Tests  of  dummy  models  proved  that  20°  and  30°  cones  could 
be  "started"  successfully  at  the  lowest  desired  stagnation  pressure 
if  the  base  diameters  of  the  models  were  decreased  to  1  inch.  Since 
the  conical  shell  of  a  30°  half -angle  cone  with  a  base  diameter  of 
only  1"  is  itself  only  1"  long,  the  initial  Type  I  design  would  have 
resulted  in  a  model  severely  affected  by  conduction  along  its  entire 
length.  Consequently,  new  20°  sharp,  20°  blunt,  30°  sharp  and  30°  blunt 
cones  were  designed  with  0.010"  skin  thickness  including  the  region  of 
the  tip.  The  plexiglass  mandrels  were  ground  to  exacting  specifications 
and  meticulous  machining  of  the  electroplated  nickel  was  required  to 
insure  uniformity  of  the  wall  thickness  both  along  the  conical  shell 
and  around  the  spherical  cap.  The  design  permitted  the  installation 
of  thermocouples  close  to  the  tip  which  were  only  slightly  affected  by 
conduction. 

No  Type  II  3°  }   5°  ,   or  15°  models  were  constructed  since  these 
models  were  long  enough  that  the  middle  and  after  portions  of  the 
conical  shells  were  virtually  unaffected  by  the  solid  tips.  The  fact 
that  identical  results  for  the  after  portions  of  the  cone  are  obtained 
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using  Type  I  and  Type  II  models  has  been  verified  experimentally  by 
Vas   .  Refer  to  Table  II  for  the  type  and  relevant  dimensions  of  the 
various  conical  models. 

D. Mounting  and  Alignment 

The  cone  models  were  supported  on  the  rear  of  the  cylindrical 
afterbody  by  means  of  a  cylindrical  holder  which  was  attached  to  a 
copper  sleeve  fitted  to  the  model.  A  water  cooled  strut,  which  could 
be  adjusted  for  pitch  and  yaw  provided  support  for  the  holder.   The 
models  were  aligned  as  closely  as  possible  to  the  tunnel  centerline, 
but  it  was  physically  impossible  to  view  the  model  by  telescope  during 
tunnel  operation  (See  Figs.  5  and  6).  The  only  check  on  the  flow 
angularity  was  provided  by  the  four  circumferentially  spaced  thermo- 
couples located  at  90°  intervals  of  arc  midway  to  the  cone  base.   If 
data  reduction  subsequent  to  the  test  revealed  equal  heating  rates, 
the  model  was  considered  to  be  properly  aligned.  If  more  than  10$ 
variation  was  observed,  the  model  was  realigned.  From  experience 
it  was  found  that  a  negative  angle  of  attack  of  ten  minutes  provided 
near  optimum  alignment. 

Ej The  Transient  Thin-Wall  Technique 

The  local  heat  transfer  coefficient  for  a  thin-skinned  body 
can  be  determined  from  a  simple  algebraic  equation  if  certain 
restrictions  are  satisfied.  For  proper  application  of  the  thin-wall 
technique  the  time  derivative  of  temperature  within  the  skin  normal 
to  the  surface  must  be  constant,  conduction  along  the  surface  of 
the  model  must  be  negligible  compared  to  the  heat  input,  and  radiation 
losses  are  required  to  be  small.  Under  such  conditions  the  following 
simple  heat  balance  for  an  element  of  skin  may  be  written: 

1  =  *-^:  =  pv0vW  (at>  II1"10 

aw  ~w  ww 
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The  term  (dv/dA)   is  determined  from  the  model  geometry  and  is 
approximately  equal  to  the  vail  thickness  for  a  cone,  while  the 
density  and  specific  heat  of  the  wall  material  are  practically 
constant  over  the  small  temperature  range  of  the  present  experiment. 
Since  the  model  walls  were  cooled  to  an  isothermal  temperature  of 
approximately  32° F  and  the  highest  temperatures  used  for  data  purposes 
was  80°F  the  assumption  of  constant  wall  properties  is  quite  sufficient 
for  the  present  experiments . 

An  important  assumption  in  the  derivation  of  the  local  heat 
balance  for  an  element  of  thin  skin  is  that  the  temperature  normal 

to  the  surface  be  uniform,  because  satisfaction  of  this  criterion 

25 
assures  the  uniformity  of  ^T/^t.  Manos  and  Taylor   point  out  that 

jL 

(ot)2>  the  distance  characteristic  of  thermal  diffusion  must  be  much 
greater  than  the  skin  thickness.     More  specifically  they  have  shown 

that  if: 

2 

''  »  1  III-ll 
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then  the  time  derivative  of  skin  temperature   in  the  direction  normal 

to  the   surface  is  constant.     Inserting  the  values  of     a  =  1.65  x  10 

o 
ft  /sec,  6   =  0.000833  ft.  and  choosing  arbitrarily     t  =  0.1  sec;  we 

obtain: 

2 
2-2p.  =  234  »  1  IH-12 

5 

which  easily  satisfies  the  inequality. 

Because  the  local  heat  input  in  the  present  experiment  is 
unsteady,   it  can  also  be  shown  that  the  Biot  number,      ho/k     must  be 
very  small  to  insure  that  the  skin  temperature  does  not  depart  signi- 
ficantly from  its  initial  value  before     dT/H     becomes  constant  in 
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26 
the  direction  normal  to  the  skin  surface.  George  and  Reinecke 

require  the  condition  that: 


ho/k  «  —  III -13 


Choosing  a  typical  value  of  h  for  a  30°  cone  test  and  inserting 
values  of  k  and  h  (h  =  2.0  x  lo"-3  BTU/ft2-Sec°R,  k  =  955  x  10 
BTU/ft-sec-°R   ,  5   =  O.OOO833  ft)  one  obtains: 


7^  =  l.Tlt  x  10"    «  0.1  III  -Ik 

k 


Consequently,  it  is  more  than  sufficient  to  establish  0.1  sec.  as 
the  time  after  which  a  uniform  dT/cH  will  be  maintained  in  the 
nickel  skin. 

The  upper  bound  on  the  time  during  which  data  can  be  recorded 

is  governed  by  the  necessity  of  limiting  diffusion  time  for  heat 

25 

tangential  to  the  model  surface.     Manos  and  Taylor       require  the 

condition  that  the  Fourier  number  based  on  a  characteristic  surface 
length  be  small  enough. to  minimize  diffusion  along  the   surface : 

2 
^    t     «    1  •  111-15 

I 

Choosing     g,   =  l/6  ft.  and     t  =  \  sec: 

2 
^-    t     =     O.306  <  1  III-16 

&d 

The  fact  that  this  inequality  is  not  sufficiently  satisfied  for  so 
short  an  upper  bound  on  testing  time  indicates  that  conduction  errors 
will  be  present.  However,  as  it  will  be  shown  in  the  section  on 
error  analysis,  the  errors  arising  from  conduction  are  not  of  significant 
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magnitude  in  the  Type  II  models.  The  conduction  problem  in  the  Type  I 
models  cannot  easily  be  calculated  due  to  the  discontinuities  in  wall 
thickness,  material,  and  metal  properties  resulting  from  the  stainless 
steel  tips. 

Several  conclusions  should  be  drawn  from  the  above  discussion. 
Moreover  it  is  obvious  that  a  well-defined  lower  bound  can  be  imposed 

on  the  testing  time  and  that  the  upper  bound  is  limited  by  conduction 

26 
errors  tangential  to  the  surface.  George  and  Reinecke   have  shown 

that  conduction  errors  may  be  expressed: 

h  at2 

Since     h     can  be  adequately  approximated  by     h  «  l/Jx     over  an  extensive 
Reynolds  and  Mach  number  range  (see  Results  and  Discussion),  the 
equation  reduces  to: 


E     =    ^    Ot  ni-18 

k        2 


where  x  is  the  distance  from  the  model  tip.  Qualitatively,  one  can 
expect  a  linear  increase  of  conduction  error  with  time,  with  the 
largest  errors  occur ing  in  the  tip  region.  In  regard  to  the  model 
design  one  observes  that  the  thin  nickel  wall  satisfies  the  major 
design  requirements.  It  should  be  noted  that  a  metal  with  lower 
thermal  diffusivity  would  satisfy  the  requirement  for  temperature 
uniformity  normal  to  the  surface  and  significantly  reduce  conduction. 
However,  materials  in  this  category  such  as  stainless  steel  are  not 
so  easily  machined,  their  thermodynamic  properties  are  difficult  to 
determine  accurately,  and  there  can  be  inhomogenities  in  the  metal. 
All  things  considered,  the  electroformed  nickel  shell  seems  to 
constitute  the  most  suitable  design  and  material. 
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p. Heat  Transfer  Apparatus  and  Data  Reduction 

The  output  of  the  copper-constantan  thermocouples  was  fed  into 
Philbrick  3&5  A  linear  operational  amplifiers,  which  were  housed  in 
a  temperature  oven  maintained  at  a  constant  100° F  to  prevent  zero  drift. 
The  gain  of  the  amplifiers  (either  25  or  50  in  the  present  work)  was 
set  high  enough  to  achieve  full  scale  deflection  of  the  60mv  galvano- 
meters in  one  second  or  less.  An  optical  arrangement  incorporated  in 
a  Midwestern  oscillograph  permitted  a  visual  representation  of  the  galva- 
nometer deflection  on  photo -sensitive  paper.  The  time  lag  of  the 
200  cps  amplifier-galvanometer  system  (order  of  50ms)  is  so  brief  that 
it  has  no  bearing  on  the  present  experiment.  A  typical  temperature  vs. 
time  trace  is  presented  in  Figure  7« 

Spot  calibrations  of  the  thermocouples  were  performed  at  room 
temperature  and  the  boiling  point  of  water.  The  output  voltages  were 

measured  with  a  hand  potentiometer  and  voltages  compared  to  National 
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Bureau  of  Standards  data  .  Within  the  three -figure  accuracy  of  the 

Rubicon  hand  potentiometer,  agreement  was  excellent,  justifying  the 
utilization  of  the  e.m.f.  vs.  temperature  relations  of  the  !NBS 
tables. 

At  the  beginning  of  each  test  the  model  was  covered  with  a 
water  cooled  copper  shield  and  maintained  at  an  isothermal  wall 
temperature  of  approximately  32° F.  The  flow  was  suddenly  established 
around  the  model  in  0.03  sec.  by  the  retraction  of  a  pneumatically 
operated  shield.  Depending  on  the  heating  rate,  from  0.5  to  1.5  sec. 
of  data  was  taken  for  record  purposes.  The  model  was  then  quickly 
cooled  to  32° F  by  admitting  cold  nitrogen  gas  into  the  conical  shell. 
Each  shell  was  vented  directly  into  the  tunnel  where  the  very  low 
ambient  pressures  effected  a  rapid  evacuation  of  the  gas  when  the 
coolant  valve  was  closed.  This  procedure  was  repeated  from  l6 
to  35  times  during  each  run  to  provide  a  large  statistical  sampling 
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for  each  thermocouple.  Every  data  point  for  an  individual  tunnel  run 
was  calculated  from  an  average  of  four  separate  thermocouple  readings 
differentiated  at  five  points.  Two  tests  at  each  stagnation  setting 
were  performed  as  a  further  indication  of  the  repeatability. 
Differentiation  of  the  data  was  accomplished  without  laborious  hand 
calculations  by  means  of  a  Model  J  OSCAR,  the  output  of  which  was 
digitized  by  an  IBM  keypunch. 

G. Error  Analysis 


Heat  transfer  data  are  subject  to  numerous  and  large  sources 
of  random  and  systematic  error.  Throughout  the  present  study  careful 
attention  has  been  paid  to  possible  sources  or  error  which  were 
prevented,  minimized  or  corrected  wherever  possible.  Unfortunately 
many  errors  arise  from  sources  beyond  the  control  of  the  investigator. 
The  major  problems  encountered  and  the  corrective  steps  taken  are 
enumerated  below. 

Systematic  Errors 
1.  Conduction 

For  a  thin -walled  heat  transfer  model,   lateral  conduction  will 
arise  not  only  from  uneven  aerodynamic  heating  but  also  from  initial 
non -uniformity  in  the  skin  temperature.      The  primary  effect   is  a 
reduction  in  measured  heating  rates  at  points  of  maximum  heating. 
Equation  III-18  demonstrates  that  conduction  effects  increase  with 
elapsed  time  and  decrease  with  distance  from  the  leading  edge . 

The  mathematical  form  of  the  conduction  term  can  be  demonstrated 
directly  from  the  simple  heat  balance  for  a  differential  element  of 
model  skin. 

dT 

p°w  6  SF  =  h  (Taw"Tw)  +    V^VT)  II1"19 
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2 
Unfortunately,  direct  evaluation  of  the  conduction  term  requires  V  T 

which  cannot  he  calculated  satisfactorily  without  closely  spaced 

thermocouples  surrounding  each  point  where  a  correction  is  to  be  made. 

Furthermore,  random  errors  in  data  reduction  render  determination  of 

the  second  derivative  most  inaccurate. 

A  simpler  method  for  estimating  conduction  effects  has  been 

28 
provided  by  Conti   whose  method  avoids  the  necessity  for  calculating 

the  Laplacian  of  temperature  by  assuming  an  aerodynamic  heating 

_i 

distribution  of  the  form  x  2  .  His  analysis  is  restricted  to  one- 
dimensional  heat  conduction,  an  initially  isothermal  surface,  constant 
adiabatic  wall  temperature,  and  negligible  radiative  heat  transfer. 
Since  each  of  these  conditions  was  fulfilled  in  the  present  experiment, 
Conti's  method  was  utilized  for  an  estimation  of  the  conduction  effects. 
Repeated  applications  of  Conti's  approximate  solution  to  typical 
data  revealed  that  most  conduction  errors  amounted  to  slightly  less 
than  1$.  For  example,  it  was  found  that  0.3  sec.  after  retraction  of 
the  boots,  the  conduction  error  had  reached  about  0.8$  at  0.6"  from  the 
model  (15°  )  tip.  In  only  the  most  limiting  cases  was  conduction  a 
significant  cause  of  error  by  an  elapsed  time  of  t  =  0.3  sec.  The 
forward  station  of  the  30°  cone,  located  0.1"  from  the  cone  vertex, 
was  the  single  most  limiting  case.  Here  the  conduction  error  reached 
12.5$  at  t  =  0.3  sec.  In  such  severe  instances  the  error  was  reduced 
to  about  5$  by  measuring  the  slope  of  the  temperature  vs .  time  trace 
during  the  initial  0.15  sec.  Cases  of  severe  conduction  could  be 
recognized  easily  by  the  deviation  of  the  trace  from  a  straight  line. 
No  actual  calculated  corrections  were  deemed  necessary  since  the 
temperature  vs.  time  data  was  differentiated  as  close  to  the  instant 
of  model  exposure  as  the  instrumentation  response  time  allowed.  Due 
to  this  procedure  transverse  conduction  losses  were  usually  less  than 
1$  but  never  more  than  about  5$. 
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An  alternate  method  of  treating  conduction  effects  has  "been 

29 

proposed  by  Thomas  and  Fitzsiramons  '    whose  method  does  not   require  an 

assumption  of  the  mathematical  form  of  the  heat  flux  distribution. 
The  correction,  which  acoounts  for  non-isothermalities  as  well  as 
uneven  aerodynamic  heating,   requires    a    least  squares  curve  fitting 
to  the  time   history  of  measured  heat  transfer  coefficients.     Although 
this  so-called  "T-F  Method"  is  more  general,  the  Conti  approach  was 
preferred  for  the  present  experiment  since  the  heat  flux  distribution 
is  known  to  closely  approximate  a  half  power  law  as  required  by  Conti' s 
analysis. 

2.  Heat  sink,  effect  of  thermocouple  wires 

Due  to  the  fact  that  the  model  skin  thickness  (0.010")  and  the 
thermocouple  wire  diameter  (0.005")  were  of  comparable  dimensions, 
it  was  conceivable  that  the  heat  conduction  along  the  wire  might 
cause  considerable  error  to  be  introduced  into  the  measurement  of 

dT  /dt    .     However,  an  estimation  of  this  error  according  to  the 

30 

method  of  McMahon   showed  the  magnitude  of  the  temperature  vs.  time 

slope  to  be  underestimated  by  2.2$  for  a  typical  heat  flux  and  t  =  1 
sec.  Although  the  error  introduced  is  shown  to  be  proportional  to 
the  thermocouple  wire  diameter  squared,  the  use  of  smaller  wire  would 
have  seriously  complicated  model  construction. 

3.  Model  skin  thickness 

Model  skin  thickness  was  carefully  controlled  in  manufacture 
and  locally  examined  to  insure  that  it  measured  0.010".  The  best 
estimate  of  average  error  in  machining  would  be  ±  .0005"  in  the 
temperature  vs.  time  history  and  therefore  ±  5$  error  in  measured 
heat  transfer  coefficients. 
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h.   Inaccuracy  in  the  knowledge  of  the  thermodynamic  properties 
of  the  electroplated  nickel 

29 

According  to  Nagel  et.al.   ,  the  specific  heat  of  nickel  has 

been  measured  to  ±  3$  for  electro formed  nickel  including  variations 

due  to  slight  differences  in  purity.  The  specific  heat  is  also  a  function 

of  temperature,  but  for  a  change  in  T   of  10°  R,  c  increases  only 

29 
0.6$   at  about  460°R. 

5.  Tunnel  axial  and  transverse  gradients 

The  reference  probe  measuring  P   was  positioned  at  the  station 
of  the  model  tip  so  that  no  correction  was  needed  for  axial  Mach  number 
gradient.  However,  corrections  were  made  for  the  transverse  tunnel 
gradient  due  to  the  fact  that  the  pitot  probe  was  5"  off  centerline. 

6.  Source  flow  effect 

Due  to  the  overwhelming  theoretical  and  practical  problems 
involved  in  the  design  of  a  contoured  nozzle  for  a  hypersonic,  high 
stagnation  temperature  wind  tunnel,  the  Princeton  nitrogen  facilities 
are  equipped  with  conical  nozzles.  As  a  result,  the  expansion  process 
continues  into  the  test  section  even  though  partial  contouring  is 
provided  by  the  tunnel  boundary  layer.   A  conical  rather  than  uniform 
free  stream  results  with  a  shock  that  is  displaced  outwards,  lower 

surface  pressures,  and  reduced  heat  transfer  and  skin  friction  coeffi- 

12 

cients   .  Although  no  data  is  available  for  the  cone,  an  estimate  of 

the  magnitude  in  the  reduction  in  heat  transfer  for  a  conical  flow  has 
been  shown  to  be  20$  for  a  wedge  at  5  inches  from  the  tip.  For  a 
conical  geometry  where  the  shock  is  closer  to  the  body  the  effect  is 
most  likely  of  at  least  the  same  order.   Fortunately  the  tunnel  boundary 
layer  provides  considerable  contouring  of  the  nozzles  used  in  the 
present  investigation.  The  flows  which  develop  more  closely  resemble 
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uniform  free  streams  rather  than  source  flows,  and  axial  Mach  number 

22 

gradients  have  "been  shown  to  be  only  0.3  per  inch   .  Although  the 

quantitative  effect  of  this  gradient  cannot  be  estimated  or  corrected 
for,  the  measurements  reported  herein  are  consistent  for  the  tunnel 
conditions. 

Random  Errors 

1.  Tunnel  effects 

The  major  uncertainty  in  test  conditions  results  from  the 

inference  of  stagnation  temperature  from  the  mass  flow  rate  (See  Sec.  IIIB) 

Due  to  scatter  in  the  flow  meter  calibrations  and  to  slight  variations 

in  the  effective  orifice  area,  the  accuracy  in  T   was  limited  to 

±  100°  K  or  ±  5$«  The  impact  of  an  error  in  T   is  compounded  by  the 

manner  in  which  heat  transfer  data  must  be  correlated.  In  the  first 

place,  T   effects  the  determination  of  Stanton  number  as  follows: 
o 

.  •  •        • 

P   =  2: ~  3" ~  3" 7—  ~  — -1 III -20 

H       p  U  c  (T  -T     )       p  U  (T  -T  )  m  3/2       t,>/S~ 

00  00   V     O      aW7  00  ooV    o      w'         p   T  •P+Vir, 

00  O         ^   ° 

Although  this  effect  is  relatively  small  (~2.25$),  a  more 
serious  inaccuracy  arises  in  knowledge  of  Re    ,  the  most  significant 

X,oo 

parameter  in  data  correlation-  Assuming  a  linear  temperature  viscosity 
relationship  for  the  low  free  stream  temperatures  encountered,  Re 
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Considering  complete  certainty  in  knowledge  of  stagnation  pressure, 
the  accuracy  of  a  free  stream  Reynolds  number  calculation  is  ±  11$  » 
by  far  the  largest  random  or  systematic  error  encountered. 
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2.  Human  error  in  data  differentiation 

The  effect  of  any  error  in  data  differentiation  is  judged  to  be 
negligible  due  to  the  large  statistical  sampling.  Each  point  represents 
the  average  of  from  12  to  l6  independent  differentiations  of  data 
obtained  in  a  time  frame  short  enough  to  preclude  the  possibility  of 
significant  variations  in  tunnel  conditions. 

Summary 

Since  errors  considered  are  independent  and  therefore  possibly 
cumulative  ,   the  accuracy  in  the  measurement  of  Stanton  number  is 
estimated  to  be  ±  15$  in  the  proximity  of  the  tip,  and  ±  11$  for  the 
downstream  stations.  The  accuracy  of  the  correlation  parameter, 
Re    ,  is  estimated  to  be  ±  11$. 
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IV 
RESULTS  AMD  DISCUSSIONS 

In  this  section,  the  sharp  cone  experimental  data  is  compared 

with  predictions  of  the  classical  Blausius -Mangier  solution,  of  the 

1  2 

viscous -invisicd  analysis  of  Probstein  including  transverse  curvature  , 

21 

of  Solomon -Loeb's  induced  pressure  and  transverse  curvature  analysis  , 

3 
and  of  Cheng's  viscous  layer  analysis  .  The  blunt  cone  results  are 

compared  with  the  local  similarity  analysis  of  Lees,  and  the  thin- 
shock-layer  theory  of  Cheng. 

A; Results  of  the  Pointed  Cone  Experiments 

The  classical  boundary  layer  predictions  for  heat  transfer  is 

given  by  Eq.  II-8  (note  that  \Jc    l  is  used  to  account  for  average 

temperature  in  the  boundary  layer).  The  expression  relates  the 

local  inviscid  Stanton  number  to  the  in  vise id  cone  Reynolds  number 

where 

^w 

St  =  IV-1 

c   p  u  (H  -H  ) 
c  c  r  w 

By  definition  the  recovery  enthalpy  H   is  used  in  the  normalization 

instead  of  the  free  stream  enthalpy  H    •  Since  radiation  losses 

O,oo 

render  the  correct  measurement  of  recovery  factor  almost  impossible 
in  hot  nitrogen  facilities,  it  is  assumed  that  r  =  \'Pr  and  that 


H  =\JPr   H  .  Inviscid  cone  values  were  obtained  as  described  in 
r   v    o 

Section  II-A.  Under  these  assumptions  the  experimental  results  are 
compared  to  the  Blausius -Mangier  theory  for  3°  >  5°  ,   15°  >  20°  and  30° 
half -angle  cones  in  Figures  8,  9  >   10;  H>  and  12>  respectively. 
(Representative  data  of  the  entire  Reynolds  and  Mach  number  range  is 
shown  in  each  figure.)  Results  for  the  3°  >  5°  >   and  15°  cones  generally 
fall  slightly  above  the  Blausius -Mangier  solution.  The  20°  cone  results 
are  scattered  about  the  predicted  value,  while  the  30°  cone  results 
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are  considerably  higher. 

In  these  same  figures  the  heat  transfer  results  are  compared  to 

the  first  order  corrections  to  the  boundary  layer  theory  due  to  the 

combined  effects  of  induced  pressure  and  transverse  curvature.   The 

results  of  the  Probstein    method  according  to  Eq.  11-13  are  shown 

for  each  cone.   The  application  of  Probstein' s  analysis  to  the 

slender  3°  and  5°  cones  results  in  a  large  overestimation  of  the 

heat  transfer,  but  the  agreement  is  fairly  good  for  the  15°  and  20° 

cones  where  the  data  generally  falls  between  Probstein 's  analysis  and 

the  boundary  layer  result .  Predicted  increases  to  the  Mangier  value 

for  the  30°  cone  are  small,  and  the  30°  results  lie  considerably  above 

the  first  order  corrected  solution  also.  The  fact  that  the  heat 

transfer  to  the  3°  and  5°  cones  is  so  highly  overestimated  may  be 

attributed  to  the  fact  that  the  transverse  curvature  parameter 

(A/r  (x))  (Eq..  11-11 )  is  so  large  in  each  of  these  cases.  The 

theoretical  analysis  requires  that  (A/r  (X))M  be  less  than  unity, 

because  the  solution  involves  series  expansions  of  a  small  parameter 

proportional  to  (A/r  (x))   .  For  the  3°  cone  a  typical  value  of 

O    m 

(A/r  (x))   is  8.9  while  for  the  5°  cone  it  may  be  3*05  •  Therefore, 
O    M 

the  expansion  is  invalid,  and  the  results  of  the  analysis  are  inapplicable 
to  the  3°  and  5°  results  of  the  present  experiment.  The  zero-order 
Blausius -Mangier  theory  gives  a  much  more  accurate  estimate  than  may 
be  obtained  by  applying  first  order  corrections  as  proposed  by  the 
Probstein  analysis.  However,  the  agreement  is  most  likely  fortuitous 
and  may  be  attributed  to  the  fact  that  the  data  for  the  slender  cones 
tends  to  fall  below  the  combined  boundary  layer  results  plus  first 
order  corrections.  At  the  lower  rarefactions  (W  <  3)  of  the  present 
study,  slender  cone  data  might  be  expected  to  begin  to  tend  toward 
the  free  molecular  value.  Evidence  of  this  phenomenon  will  be 
discussed  in  greater  detail  below. 
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The  heat  transfer  data  are  correlated  in  terms  of  the  parameters 

3  8 

proposed  by  Cheng  and  modified  by  Waldron  for  the  viscous  layer 

regime  in  Figs.  13,  Ik,   15,  l6  and  IT: 
C  Re 

H,oo    .       TTT/    v  X.oo        1 

*\   '  (1  -  W/Tn)  VS.  w  =  — i —  IV-2 

sine  v        °J  2  7 c  cose 

c  M     *    c 

00 

Cheng  originally  suggested  that  C   /sinec  be  used  in  the  correlation, 

H,°°  b 

but  for  simplicity     e_     has  been  replaced  by     fl      .     To  account  for  the 
variation  in  wall  and  stagnation  conditions  the  ordinates  have  been 
modified  by  the  factor  (1  -  T^/Tq)  as  suggested  by  Waldron   .     The 
constant  used  in  the  linear  viscosity -temperature  relationship  is 
defined  such  that: 


T       =     T     [1  +  3(T  /T  )]/6  IV-3 

*  o  wo 

5  3 

as  proposed  by  Cheng  .  In  his  later  viscous  layer  analysis  he 

suggested  that  the  reference  temperature  for  conical  flows  be  taken  as: 

T   =  ^v/p  +  T  (-  -  -   cos  e  )  iv-k 

*  ^    0  2   3        cy 

As  e   decreases  the  latter  reference  temperature  reduces  to  the  former, 

which  can  be  shown  to  be  the  average  temperature  across  a  locally 

similar,  hypersonic  boundary  layer  on  a  flat  plate.  No  significant 

difference  arises  in  the  two  reference  temperatures  until  9     =  20°  . 

c 

For  the  20°    and  30°    cone  of  the  present  study     c       would  be  reduced  by 
5$  and  10J&  respectively  had  Eq..   IV-k  been  used.     For  the  20°    and  30° 
cone  the  data  is  more  effectively  correlated  using  Eq..   IV-3,   hence  this 
equation  is  used  throughout. 

Even  though  Cheng's  heat  transfer  parameters  were  proposed  for 
the  viscous -layer  regime  for    W  <  10,  they  have  proven  valid  for 
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correlating  heat  transfer  data  over  a  wide  range  of  heating  rates, 

cone  angles,  and  Reynolds  numbers.  However,  in  the  present  results 

a  slight  dependence  of  C   /sin0   can  be  observed.  For  a  given 

H/»     c 

value  of  W  a  general  increase  in  C   /sine   occurs  as  cone  angle 

H,°°     c 

decreases.  This  increase  may  be  observed  consistently  except  in  the 
case  of  the  low  Reynolds  number,  30°  results  which  fall  higher  than 
the  corresponding  20°  results.  A  viscosity  relationship  based  on 
T*  as  defined  by  Eq..  IV-^-  would  not  affect  this  trend  for  the 
smaller  cone  angles  where  the  neglected  transverse  curvature  effects 
are  much  more  significant. 

A  modification  of  the  boundary  layer  theory  to  include  the 

first  order  effects  of  induced  pressure  and  transverse  curvature  has 

21 
been  proposed  by  Solomon-Loeb   .  Results  of  their  Pohlhaussen  solution 

are  presented  with  the  data  for  each  cone  in  Figure  13-17.  It  is 
important  to  note  that  the  Solomon-Loeb  solution  is  not  limited  to 
small  values  of  the  transverse  curvature  parameter  (A/r  (x))  as  was 
Probstein's  analysis.  In  fact  the  former  authors  give  a  limiting 
solution  for  infinite  transverse  curvature  effects  (Eq.  II-16)  which 
is  of  course  valid  for  very  slender  cones  at  hypersonic  Mach  numbers. 
In  direct  contrast  the  validity  of  Probstein's  results  for  slender 
cones  would  be  limited  to  supersonic  Mach  numbers  for  the  present,  low 
Reynolds  number  range. 

Agreement  of  the  heat  transfer  results  and  the  Pohlhaussen 
solution  are  excellent  for  the  15°  ,   20°  ,  and  30°  cases  throughout  the 
Reynolds  number  range.  For  the  slender  3°  and  5°  cones  the  results 
are  good  for  W  >  100  .  The  range  of  W  covered  in  the  present 
study  is  particularly  interesting  because  the  first  evidence  of  the 
departure  of  the  results  from  the  boundary  layer  theory  may  be 
observed.  The  3°  and  5°  results  appear  to  depart  from  the  Solomon- 
Loeb  solution  at  W  ~  60-90  and  tend  toward  the  viscous  layer 
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3 
analysis  of  Cheng  .  These  results  are  consistent  with  the  observation 

8 
of  Waldron  who  commented  on  a  similar  departure  of  the  results  of 

T 

Wilicingson  and  Harrington  in  the  neighborhood  of  W  =  100  .  As  the 

boundary  layer  limit  for  his  observation,  Waldron  employed  Probstein's 
first  order  transverse  curvature  correction  (Eq.  11-12).  This 
procedure  is  incorrect  for  two  reasons:   In  the  first  place  Probstein's 
analysis  is  invalid  for  (A/r  (*))„  »  1,  and  secondly  Waldron  failed 
to  consider  the  induced  pressure  correction  to  the  Mangier  solution 
(Eq..  11-10).  The  6.3°  results  of  Reference  7  actually  agree  well 
with  the  present  5°  results  as  well  as  to  the  proper  boundary  layer 
limit.  Departure  from  the  Solomon-Loeb  solution  of  the  shock  tunnel 
data  occurs  at  W  =  90  . 

For  the  larger  half -angled  cones  departure  from  the  boundary 
layer  result  does  not  occur  until  the  flow  is  more  rarefied  at  lower 
values  of  W  .  The  15°  cone  results  follow  the  Solomon-Loeb  solution 
at  least  as  far  as  W  =  15  >  while  evidence  of  departure  for  the  20° 
and  30°  cones  can  be  first  observed  at  W  ~  3  • 

Selection  of  a  particular  limiting  case  for  the  Pohlhaussen 
solution  was  made  with  several  considerations  in  mind.  The  transverse 
curvature  parameter  X  can  be  related  to  W  by  a  constant  for  a 
given  cone  angle.  This  constant  is  large  enough  so  that  T  is  not 
sufficiently  small  to  apply  the  limiting  solution  for  X  -  0  except 
for  the  3°  cone.  For  cones  of  larger  half  angles  X  remained  close 
enough  to  unity  to  cause  Eq..  II-16  to  blow  up  because  of  the 
denominator  X  log  =  .  Therefore,  the  limiting  solution  for  small 
transverse  curvature,  Eq.  11-17,  was  utilized  for  comparison  to  the 
5°  ,  15°  ,  20°  ,  and  30°  cones. 

B_. Results  of  Blunt -Nose  Cone  Experiments 

In  an  effort  to  determine  the  validity  of  several  theoretical 
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treatments  of  the  blunted  cone,  experimental  results  were  obtained 
on  spherically-capped  right  circular  cones.  A  wide  range  of  length- 
to-nose  radii  is  investigated  for  each  of  the  5°  >   15°  ,   20°  ,   and  30° 
cones. 

The  experimental  results  were  compared  to  the  theoretical 

k 
calculation  of  Lees  whose  result  for  the  heat  transfer  distribution 

on  the  conical  surface  of  a  blunt  cone  may  be  written: 

6  x/r 

p-T"  =  A(0  ) r-r-7^  IV-5 

Cqv>o    c  tB(ec)  +  (*/r0)3]1/2 

for 

—  >  cote  iv-6 

r    -     c 
o 

x  ,   the  distance  along  the  surface  measured  from  the  virtual  tip  is 
related  to  the  actual  distance  from  the  stagnation  point  by 


x  =  S  +  r  [cote  -  (n/2  -  e  )] 
o     c   *  c    c' 


The  function  A(e  )  and  B(e  )  depend  on  7 ,  M   ,  and  G      (Ref.  h) . 

x  C  C  ■  CO  '         C 

The  local  heat  flux  must  be  normalized  by  the  stagnation  point  heat 
transfer  which  was  not  measured  in  the  present  experiment.  Consequently, 
(a  )   was  calculated  from  the  Lees  expression: 

,     s         |  Pr'2/3  yrC^To  yftT   hse  (Km,,  7  ,  y  ) 

(     )    = a_ IV.T 

Results  for  the   5°  >   15°  >   20°    and  30°    cones  are  presented  in 
Figures  18-21.     The  several  test  section  conditions  represented  cover 
the  entire  Mach  and  Reynolds  ranges.     The   5°    results  show  neither 
qualitative  nor  quantitative  agreement  with  the  simplified  Lees 
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expression.  Data  points  spanning  the  entire  Reynolds  number  range  fall 
within  ±  8$  ,  it  is  clear  that  the  Lees  formulation  does  not  apply- 
to  very  small  half -angle  cones  of  the  present  bluntness.  Ratios  of 
%  '(%r)       ^or  the  15°  ,  20°  and  30°  cones  correlate  well  over  the 
range  of  conditions,  and  at  least  qualitative  agreement  with  the 
theory  is  exhibited  by  the  slopes  of  the  data.  However,  values  of 
a  /(a  )   fall  significantly  above  the  theoretical  distribution.  For 
the  15°  cone  this  discrepancy  ranges  from  25$  near  the  leading  edge 
to  45$  near  the  trailing  edge,  22-g-  nose  radii  downstream.  The  predicted 
heat  transfer  is  underestimated  by  approximately  35$  to  55$  for  the  20° 
cone,  and  by  about  24$  to  35$  for  the  30°  cone.  Due  to  the  consistent 
underestimation  of  the  ratio  q  /(q  )   ,  one  must  examine  the  possibility 
of  an  underestimation  of  the  stagnation  point  heat  transfer.  By 

comparison  of  the  Lees  result  to  the  experimentally  verified  solution 

31 
of  Fay  and  Riddell  ,  one  finds  agreement  accurate  to  less  than  3$« 

It  is  therefore  unlikely  that  the  calculated  values  of  (a  )   were 

seriously  in  error. 

The  experimental  blunt  cone  results  were  compared  to  the  zero 

5 
order  predictions  of  Cheng  who  proposed  the  correlation  parameter 

m  2(€k)lAc       e2  s 

oo  H,oo  C 

VS. 


6  2T  Vek  t 

c   t 

The  domain  of  validity  of  Cheng's  formulation  is  characterized  by 

several  restrictions  which  govern  the  degree  of  accuracy  of  the 

analysis   in  the   inviscid-flow  region  as  well  as  the  boundary  layer, 

These  conditions,  which  were  discussed  above,  may  be  summarized 

and  examined  as  follows: 
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Theoretical  Assumption 


1)     7 


7-1 


2)  c    i  '    -/7+1  «  1 

3)  (1/2)0     2^"1)^  «  i 


2       2 
4)     M       0       »  1 

oo  C 


Experimental  Condition 

7    -  1.1* 

e  =  1/6 

.124  <(l/2)0c  2^"1)/7  <   .345 

2        2 
M       9        =2.2     for     5 
00         c 

19.8  for  15° 
35.2  for  20° 
79 -2  for  30° 


The  several  restrictions  listed  are  mutually  exclusive  in  that 
they  require  a  large  flow  deflection  angle  to  satisfy  (k)   and  at  the 
same  time  a  small  deflection  angle  to  satisfy  (3)-   The  latter  two 

conditions  are  more  nearly  satisfied  by  the  15°  cone.  The  limit  of 

2   2 
validity  of  the  strong  shock  condition  (M  9       »  l)  is  approached 

00    c 

by  the  5  cone  while  the  limit  of  Cheng's  small  perturbation  velocity 

2(  y  -1  /y  ) 

requirement  (l/2  Q  «  l)   is  approached  with  the  30°    cone. 

Results  of  the  comparison  of  the  data  to  Cheng's  zero-order 
shock  layer  theory  are   shown  in  Figures  22-25.      The   nose  drag 
coefficient,   k   ,  used  in  the    correlation  was  calculated  according  to 
a  modified  Newtonian  pressure  distribution  for  a  spherically  blunted 
nose  (See  e.g.  Ref.   32).     The  numerical  agreement  of  the  data  with 
the  theory  is  generally  poor,  and  the   power  law  slope  predicted  is 
incorrect.      However,   the  correlation  parameters  themselves   seem  to 
be  quite  good  since  the  data  for  the  full  range  on  cone  angles  and 
Reynolds  numbers  follow  the  same  power  law  curve  (See  Fig.   26).     The 
poor  agreement  of  theory  and  experiment  is  most  likely  due  to  a  break- 
down of  the  assumptions  of  the  analysis  and  the   zero-order  nature  of 
the  solution.     Where  the   small  perturbation  velocity  assumption  breaks 
down  for  the  larger  cones,   the  heat  transfer  is  underestimated;  and 
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where  the  shock  weakens  for  cones  of  small  half -angles,  the  heat 
transfer  is  underestimated. 

The  predicted  oscillatory  decay  is  evidently  associated  with  the 
zero-order  nature  of  the  solution  and  does  not  have  the  physical 
significance  of  a  skipping  and  reattaching  shock  layer  as  Cheng 
suggested.  No  oscillations  in  the  data  are  observed  in  the  results 
for  any  cone  angle . 

C_. Comparison  of  Sharp  and  Blunt  Cone  Results 

Due  to  the  high  heat  transfer  rates  encountered  in  the  stagnation 

regions  of  actual  flying  vehicles,  at  least  some  nose -"blunting  is  required 

to  reduce  the  heating  problem  to  manageable  proportions.  For  practical 

engineering  applications  it  is  of  interest  to  know  the  convective  heat 

transfer  rates  at  equal  distances  from  the  virtual  tip  of  a  sharp  or 

blunt  cone  of  a  given  half -angle.   In  Figures  27  ond  28  the  ratios  of 

the  blunt  cone  heat  transfer  parameters,  C  ,  ,  to  the  sharp  cone  heat 

H,b 

transfer  parameters,  C    ,   are  presented  for  the  full  Reynolds  number 
range  for  the  5°  and  15°  cones.  Since  these  cones  were  of  the  inter- 
changeable tip  variety,  it  was  possible  to  change  tips  without  affecting 
the  alignment  of  these  models.  Ratios  of  these  parameters  are  not 
presented  for  the  Type  II,  20°  and  30°  cones  for  several  reasons: 
1)  Installation  of  separate  sharp  and  blunt  models  made  it  impossible 
to  ensure  exactly  the  same  alignment  or  skin  thicknesses,  and  2)  Due 
to  calibration  problems  the  30°  cones  were  tested  at  Reynolds  numbers 
which  differed  by  more  than  10$  between  the  sharp  and  blunt  cases. 

These  uncertainties  and  the  differences  in  test  conditions  completely 

k 
obscure  effects  due  to  bluntness  alone.  However,  Lees  has  shown  that 

variations  between  heating  rates  of  sharp  and  blunt  cones  of  large 

half -angles  are  small. 

For  the  5°  cone  the  reduction  in  heat  transfer  due  to  bluntness 
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is  as  large  as  20$.  For  x  =  2"  the  blunt  cone  heating  rates  on  the 
5°  cone  are  considerably  larger  than  the  sharp  cone  rates.   This  result 
is  consistent  with  the  data  of  Feik  ,   who  tested  5   cones  of  the  same 
bluntnesses  under  similar  testing  conditions.  At  locations  in  the 
vicinity  of  x/d  -   13  the  blunt  cone  heating  rates  fall  about  20$ 
below  the  sharp  cone  results.  For  the  lower  Reynolds  number  cases 
the  bluntness  effects  disappear  in  the  neighborhood  of  x/d  =  2k   , 
while  the  blunt-cone  entropy  layer  for  the  higher  Reynolds  number 
case  still  has  considerable  influence  beyond  that  point. 

For  the  15°  cone  the  influence  of  bluntness  is  greatest  at 
the  first  data  point,  x/d  =  k    .  Here  the  blunt  cone  results  fall 
some  20$  below  the  sharp  cone  heating  rates.  The  influence  of  the 
blunted  tip  gradually  decreases  as  one  moves  further  downstream  on 
the  cone.  At  x/d  values  of  about  10,  no  further  bluntness  effects 
are  detected  in  either  of  the  higher  or  lower  Reynolds  number  cases. 
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V 
CONCLUSIONS 

The  most   important  results  of  the  present   study  are: 

(1)  The  classical  boundary-layer  solution  provides  a  good 
estimate  of  local  heat  transfer  distributions  for  the  Reynolds  and 

Mach  number  ranges  investigated  when  an  appropriate  reference  temperature 

(T   )   is  chosen  for  the   linear  temperature -viscosity  relationship. 

1  2 

(2)  The  large  rise   in  heat  transfer  predicted  by  Probstein  ' 

in  his  transverse  curvature  and  boundary-layer  displacement  analyses 

3 
is  not  observed  experimentally  even  for  the  3°    cone.     However,  Cheng's 

viscous-layer  solution  for    W  <  10     is  quite  accurate   in  this  portion 

of  the  transition  regime  between  continuum  and  free  molecular  flow. 

Although  the  continuity,  momentum  and  energy  equations  solved  by 

Probstein  and  Cheng  are  basically  the  same,   a  substantial  difference 

is  evident  in  the  numerical  results  due  to  the  fact  that  Probstein's 

outer  boundary  conditions  are  given  by  the  external  inviscid  flow 

whereas  Cheng's  outer  boundary  conditions  are  given  at  the  shock. 

Since  the  present  results  support  Cheng's  predictions,   the   conclusion 

that  transport  at  the  shock  is  important   in  reducing  heat  transfer 

below  the   classical -boundary  layer  limit  at  low  Reynolds  numbers  is 

reinforced. 

,    x  21 

(3)  The  limiting  solution  of  Solomon-Loeb     ,  which  includes 

transverse  curvature,   is  valid  over  the  entire  cone  angle  range 

(except  3°  ?   low  Reynolds  number)  for    W  <  10    . 

3 
(k)     Cheng's     parameter,  W   ,    is  valid  for  correlating  heat 

transfer  data  over  the  entire  Reynolds  number  and  cone-angle  range. 

The  experimental  data  expressed  as  (C  /sin  0    )(1-T  /t   )  exhibit  only 

Hcwo 

a  small  dependence  on  cone  angle.     The  increase  in  heat  transfer  for 
small  cone  angles  is  attributed  to  transverse  curvature;  however,   the 
rise  is  much  smaller  than  that  predicted  by  Probstein. 
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(5)  For  engineering  design  purposes  the  theoretical  analyses 

3  21 

of  hoth  Cheng  and  Solomon -Loeb   are  completely  adequate.  For  the 

higher  Reynolds  numbers,  W  <  90  }   "the  classical  boundary  layer 

solution  is  quite  satisfactory. 

k 

(6)  The  predictions  of  Lees'   simplified  blunt  cone  method 

fall  below  the  experimental  data  for  all  cone  angle  tested. 

(7)  Cheng's  bluntness  parameter  correlates  data  for  the  full 
Reynolds  and  cone-angle  range,  although  the  oscillations  predicted  by 
his  thin  shock  layer  theory  are  not  observed. 

(8)  The  smaller  the  cone  angle  the  more  effective  is  nose 
blunting  in  reducing  heat  transfer  below  the  equivalent  sharp  cone 
results.  The  5°  blunt  cone  results  exhibit  a  20$  reduction  in  heat 
transfer  as  far  as  twenty  nose  diameters  downstream.  The  15°  results 
show  a  20$  reduction  at  four  diameters  while  no  reduction  at  all  is 
evident  twelve  diameters  downstream.  The  20°  and  30°  cone  data 
exhibited  heat  transfer  results  which  were  insensitive  to  nose 
blunting  less  than  one  nose  diameter  downstream. 

The  following  recommendations  are  provided  for  extensions  of 
the  present  study: 

(1)  The  uncertainty  in  wind  tunnel  conditions  could  be  sub- 
stantially reduced  if  the  graphite  heating  elements  were  pre -treated 
to  eliminate  carbon  flaking  and  subsequent  deposits  in  the  nozzle 
throat. 

(2)  Heat  transfer  is  sensitive  to  stagnation  conditions. 
Although  stagnation  pressure  can  be  determined  precisely  with  Heise 
gages,  stagnation  temperature  is  difficult  to  determine  accurately. 
Meticulous  care  should  be  taken  in  the  calibration  of  nitrogen 
facilities  for  heat  transfer  experiments. 

(3)  Further  heat  transfer  investigations  are  needed  at  lower 
Reynolds  numbers  than  could  be  obtained  in  the  present  facilities 
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in  order  that  the  transition  regime  "between  continuum  and  free- 
molecular  flow  can  be  more  fully  understood. 
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TABLE  II 


Half  Angle  Based  Length 

9  Type       R      inches      (inches  from  virtual  apex) 


002       .750  7. IT 

1.50  8.59 


.002 
.1125 


10*  I  '™  1.50  4.32 

15°  I  [™%  1.50  2.90 

20°  II  .002  1.00  1.46 

20°  II  .1125  1.00  1.46 

30°  II  .002  1.00  1" 

30°  II  .1125  1.00  1" 
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Figure  T«  Typical  Oscillograph  Traces  of  Voltage  vs.  Time 
for  the  20  Spherically  Blunted  Cone. 
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3°    Sharp  Cone  Stanton  Number  vs.  Reynolds  Number 
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5°  Sharp  Cone  Stanton  Number  vs.  Reynolds  Number 
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Figure  10 


15°  Sharp  Cone  Stanton  Number  vs.  Reynolds  Number 
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Figure  27 
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Figure  28 


Ratio  of  Blunt  Cone  to  Sharp  Cone  Heat  Transfer 
Parameters  for  the  15°    Cone  at  Equal  Distances 
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